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ABSTRACT 



Context. In close binary systems, the axial rotation and the mutual tidal forces of the component stars deform each 
other and destroy their spherical symmetry by means of the respective disturbing potentials. 

Aims. We present new models for low-mass, pre-main sequence stars that include the combined distortion effects of tidal 
and rotational forces on the equilibrium configuration of stars. Using our theoretical results, we aim at investigating 
(— I ' the effects of interaction between tides and rotation on the stellar structure and evolution. 

Methods. The Kippenhahn & Thomas (1970) approximation, along with the Clairaut-Legendre expansion for the grav- 
I ■ itational potential of a self-gravitating body, is used to take the effects of tidal and rotational distortions on the stellar 

O ' configuration into account. 

^ ' Results. We obtained values of internal structure constants for low-mass, pre-main sequence stars from stellar evolu- 

C/3 tionary models that consider the combined effects of rotation and tidal forces due to a companion star. We also derived 

a new expression for the rotational inertia of a tidally and rotationally distorted star. Our values corresponding to 
standard models (with no distortions) are compatible with those available in literature. Our distorted models were 
successfully used to analyze the eclipsing binary system EK Cep, reproducing the stellar radii, effective temperature 
ratio, lithium depletion, rotational velocities, and the apsidal motion rate in the age interval of 15.5-16.7 Myr. 
Conclusions. In the low-mass range, the assumption that harmonics greater than j=2 can be neglected seems not to 
be fully justified, although it is widely used when analyzing the apsidal motion of binary systems. The non-standard 
evolutionary tracks are cooler than the standard ones, mainly for low-mass stars. Distorted models predict more mass- 
concentrated stars at the zero-age main-sequence than standard models. 



> 

in 
cn 
in 

(N 

^' 
O 

a^ 
o 

> 

>< 



Key words. Stars: evolution - Stars: interiors - Stars: rotation - Stars: pre-main sequence - Stars: binaries: close - 
Stars: low-mass, brown dwarfs 

1. Introduction to the Kcplcrian angular velocity, ujk, of the orbital motion 

around a common center of mass, so that 



A binary system consists of two stars that rotate around 

■ their own axes and, at the same time, revolve around the i — 

5^ [ center of mass of the system. Sometimes, the intrinsic rota- fl ^ ujj^ = J G—^ —. (1) 

tion axis and the orbital one are aligned since the beginning ^ 
of the formation process. Usually, the orbit begins with a 

considerable eccentricity and the component stars are not However, occasionally il is much larger than ujk - the 

synchronized with the orbital angular velocity. However, sense of rotation bemg direct m every known case. In sys- 

due to the inertial forces that take place, the system tends ^^ms exhibiting circular orbits, synchronism between rota- 

to align its axes, to synchronize the rotational angular ve- ^^O" ^nd revolution may usually (though not always) be 

locity of the components with the orbital angular velocity at expected to exist, while components describing eccentric 

the periastron passage and, finally, to circularize the orbit, ^^'bits, though m general rotating faster than their mean 

Extensive spectroscopic evidence reveals that the compo- o^^ital angular velocity, can still be synchronized with the 

nents of close binary (in general, non-eccentric) systems do orbital motion at the periastron. 

rotate with an angular velocity, O, which is generally equal When both components arc^ on the main sequence, the 

most massive one is also larger in radius and hotter, and is 
called primary, both in photometric (the star eclipsed in the 



Send offprint requests to: N.R.Landin deepest minimum of the light curve) and in spectroscopic 

* The complete version of Table[T] is only avail- studies. In other evolutionary stages, however, it is a com- 

able in electronic form at the CDS via anony- mon situation that the most massive component is not the 

mous ftp to cdsarc.u-strasbg.fr (130.79.128 .5) or via larger or the one of higher effective temperature. It contin- 



http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-|-A/[ ues being designed as "primary" in spectroscopic studies. 
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but referred to as "secondary" in light curve analysis, due 
to the lower effective temperature. 

One important aspect of the evolution of close binaries 
is the dynamical evolution due to tidal interaction, which 
is reflected in the rotation of the stars and in the eccentric- 
ity of their orbits. Tidal deformation due to the companion 
would be symmetric about the line joining their centers, if 
there were no dissipation of kinetic energy into heat. It is 
this dissipation that induces a phase shift in the tidal bulge, 
and the tilted mass distribution, then, exerts a torque on 
the star, leading to an exchange of angular momentum be- 
tween its spin and the orbital motion. Theory distinguishes 
two component s in the tide , namely, equilibrium tide and 
dynamical tide (|Zahnl . Il989t ): 

— Equilibrium tide is the hydrostatic adjustment of the 
structure of the star to the perturbing force exerted by 
the companion. The dissipation mechanism acting on 
this tide is the interacti on between the convective mo- 
tions and the tidal flow (jZahnl . Il966f ) . 

— Dynamical tide is the dynamical response to the tidal 
force exerted by the companion; it takes into account 
the elastic properties of the star, and the possibilities of 
resonances with its free modes of oscillation. The dissi- 
pation mechanism acting on this tide is the departure 
from adiabaticity of the forced oscillation, due to the 
radiative damping (jZahnl . Il975[ ) . 



To describe the tidal process in massive main-sequence 
stars, which have convective cores and radiative envelopes, 
it is necessary to use a theory that account s for dynam- 
ical effects which aris e due t o tidal forces (|Zahnl. 11977 ; 
Savoniie fc Papaloizoul Il983t) . IWitte fc Savoniid (|l999al lb 



200 it) published extensive studies of the role played by close 
resonances with eigenmodes of these early type s tars dur- 
ing the orbital decay (see also lWillems et an . l2003[ ). On the 
other hand, for low-mass main-sequence stars and giants, 
which have extended convective envelopes, retardation of 
the equilibrium tide due to the viscosity of turbulent eddies 
in the envel ope is usually assumed to be the cause of the 
tidal torque. [Savoniie fc Wittd (|2002a[ ) calculated the tidal 
interaction of a uniformly rotating 1 Mq star with a orbiting 
companion at various phases of evolution from the zero-age 
main-sequence (ZAMS) to core hydrogen exhaustion. Their 
results indicate that effects related to stellar rotation can 
considerably enhance the speed of tidal evolution in low- 
mass b inary systems. In another paper, IWitte fc Savoniid 
(|2002bD showed that energy dissipation through resonant 
dynamic tides may dominate convective damping of equi- 
librium tides in solar-type stars. 

In standard models, the stars are assumed to be spher- 
ically symmetric. However, the spherical symmetry will be 
destroyed if a disturbing potential exists, as it happens for 
rotating stars either isolated or in binary (multiple) sys- 
tems, or due to tidal forces (gravitational influence of a 
companion), present in binary systems. 

In the case of a rotating star in a binary system, 
both rotational and tidal forces distort its shape from the 
spherical symmetry. The analytic determination of these 
combined effects is complex and approximate methods 
have been used in the literature, with one of the distort- 
ing forces (g enerally rotat i on) be i ng analyzed in approx - 
imate ways () Mohan et al.l . Il990f ). IChandrasekhail (I1933D 
developed t he theory of distorted polytropes and iKopall 
()l972l IT974I ) developed the concept of Roche equipoten- 



tials and coordinates to study t he combined effects of 
tidal forces and rotation on stars. iKippenhahn fc ThomasI 
(jl97(j ) (KT70) devised a method for introducing the ef- 
fects of rotation in existing one-dimensional evolution- 
ary codes tha t became wid el y ad o pted in the litera- 
ture ( see e.g. Endal fc Sofial. 19761: Pinsonneault et al. 



119901; iFliegcner fc Langcii 119951: iChabover et all 11995 
Mevnct fc Maeder. . ,199'fl : iMendes et all . Il999at IClaret 
1999D . In the case of massive stars, the effects of rotation on 



their struc ture and evolution were studied by various au- 
thors (e.g. iMevnet fc Maedeij|l997l. l200Gl: iMaeder fc ZahnI 



119981: iHeger et al.l l2000t iMaeder fc Mevnetl l2000l l2001li 

some of which present also grids of stellar models for ro 
tating, massive stars. 

In the low-mass r ange , Pinsonneault et al. _( 1990[ ) 
iMartm fc Claretl (|1996[ ) and (Mendes et al.l (|1999ah pre 



sented exploratory studies on the influence of rotation on 
Li depletion i n pre -main sequence (MS) evolution, and 
iLandin et al.l ()2006f ) provided grids of non-gray rotating 
pr e-MS models f r om 0. 085 Mfn to 3.8 Mg. 

iMohan et al.l ([19901 ) presented a method for calculat- 



ing the equilibrium structure of a rotationally and tidally 
distorted primary component in a synchronously rotating 
binary system. Their method assumes a Roche model po- 
tential and uses the KT70 method for computing the cor- 
rections due to rotation and tides. Their results for stellar 
models of 10, 5, and 2.5 Mq, for which the mass ratio q of 
the secondary component to the primary component was 
set to 0.1, indicate that rotational effects are more impor- 
tant than those produced by tidal distortions. 

In this paper, we p resent a new version of the ATON code 
(j Ventura et al.l . Il998t ) that treats the combined effects of 
tidal and rotational distortions on a star. This allows us, for 
example, to obtain the apsidal motion constants of the pri- 
mary component of a binary system that rotates and suffers 
the tidal effects of the secondary component. Instead, how- 
ever, of using a Roche-type potential, which considers both 
binary components as point masses, we adopt the more pre- 
cise technique based on the Clairaut-Legendre ex pansion 
for th e potential energy of a self-gravitating body ijKopall . 
|1959[ ). The paper is organized as follows: in Sect. ^ we 
give a brief account of theoretical apsidal motion calcula- 
tions in the literature. Sect. ^ presents the KT70 method, 
proposed for determining the equilibrium structures of ro- 
tationally and tidally distorted stellar models, in which the 
non-spherical stellar equations can be easily obtained from 
the spherical ones. The technique for introducing the com- 
bined tidal and rotational effects is described in Sect. in 
which wc present, also, new calculations of internal struc- 
ture constants extended to the pre-MS phase. The results 
are presented in the Sect. ([5]). Discussion and comparisons 
with observed apsidal motion rates are given in Sect. ©. 

2. Apsidal motion and internal structure constants 

The internal structure constants k2, and k^, also known 
as apsidal motion constants, are important in stellar astro- 
physics. They are mass concentration parameters that de- 
pend on the mass distribution throughout the star. There 
is, however, a direct relation between the gravitational field 
of a non-spherical body and the internal density concentra- 
tion in that body (jSahade fc Wo od. 1978[) . 

From a theoretical point of view, the values of kj 
(j=2,3,4) depend on the model used. For the Roche model. 
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in which the whole stellar mass is concentrated at its cen- 
ter, the kj values are all equal to zero, while for a homo- 
geneous model ^2=3/4, fc3=:3/8 and k4=l/A. The values of 
the internal structure constants are essential to compute the 
theoretical apsidal motion rates in close binaries, and the 
comparison with the observations constitutes an important 
test for evolutionary models. The most centrally concen- 
trated stars have the lowest values of kj and the longest 
values of apsidal periods (Eq.[6]). 

It can be shown that the theoretical apsidal rate uj, in 
radians per cycle, is given in terms of the in t ernal struc- 
ture constants by ijMartvnovl . Il973l : iHeileserJ . Il987l . here- 
after H87) 



2^ 



EE 

i=l 3=2 



Cjikji , 



(2) 



where i denotes the component star (l=primary, 2=sec- 
ondary) and j the harmonic order. Generally, terms of or- 
der higher than j=2 are very small so that only values of 
k2i enter in Eq. ([2]). 

The internal structure constants are important in other 
astrophysical aspects, since synchronization and circular- 
izatio n time scales in close binaries depend on ^2 (jZahnl . 
Il977t) . Other applications are in the computation of rota- 
tional angular momenta (Sect. H^ . where gyration radii 
(defined in Eq. [58|) can be ex pressed as a lin ear function of 
the apsidal motion constants (jUrechel . ll976l ). and in the de- 
termination of the effect of binarity in the geometry of the 
stellar surfac es due to rotation and tides (.Ruciiiskii . Il969t 
lKopa]| . fl97l . 

The first analytical expression for the apsidal motion 
period in close binaries in terms of stellar masses, rela- 
tive radii and internal structur e cons tants of the compo- 
ne nt stars was giv e n by Russelll fl928r) an d later improved 
by ICowHng] (|1938[ ) . IChandrasekhail ( 1933t) used polytropic 
models to predict internal structure constants for main- 
sequence stars. At that time, the large uncertainties of the 
observational data, as well as the use of polytropic mod- 
els with an arbitrary index n, were responsible for the ap- 
parently good agreement between observed and predicted 
values of log ^2 . 

By using more realistic stellar models, different au- 
thors derived more elaborated expressions for the apsidal 
motion period, separating rotational and tidal contribu- 
tions to the total apsidal motion rate. The apsid al mo- 
tion t e st was also applied to polytropic models by ISterni 
(|1939f) . lBrooker fc Olid (|1955f) and, l ater, to early t h eoret - 
ical stellar models a t the ZAMS, bv ISchwarzschildl (I1958D 
and Kushawal (jl957f ) , both using the old Keller fc Meverotg 
()1955h opacities. At that time, the theoretical values of ^2 
were systematically larger than those obtained by observa- 
tions, which led to the interpretation that real stars were 
more centrally condensed than predicted by models. This 
discrepan cy persisted during several deca des; more recently, 
howev er, IClaret fc Wi llcms {2QM) and iWillems fc ClarelJ 
()2003f ) have found no systematic effects in the sense that 
models are less mass concentrated than real stars, for bi- 
naries with accurately known absolute dimensions. iJeffervl 
()1984f ) and H87 computed internal structure const ants for 
stars within the main sequence. The former used ICarsonI 
(Il976l) opacitie s , whi le the latter used opacity tables by 
iCox fc StewartI (Il969h . 



In more recent years, the most extensive series of the- 
oretical works in the literature regarding apsidal motion 
constants is that one from A. Claret and collaborators, 
some o f which are briefly described here. lClaret fc Gimened 
(|l989al . hereafter CG89a) presented a detailed grid of evolu- 
tionary stellar models during the hydro gen burning phases, 
i nclud ing apsidal motion constants. IClaret fc Gimenea 
(Il99ll ) studied the effect of the core overshooting and mass 
loss on the internal density concentratio n of main-sequence 
stars. With more updated input physics. fClaret fc GimenezI 
(|1992| . hereafter CG92) computed stellar models together 
with internal together with internal structure constants. 
Those models showed a general tendency to be cooler and 
more centrally concentrated in mass than their previous 
computation s for the typical ma s ses w h ere ap sidal motion 
is observed. IClaret fc GimenezI (| 19921 \l99m and IClaretl 
(|1995D were able to reduce the discrepancies between 
theoretical and observed values of ^2 at acceptable levels 
for sys t ems w hose relativistic contributions were small. 
IClaretl (|l999t ) took into account the effect of rotation 
on the internal structure of stars and found that it 
strongly depends on th e distortion of the configuration. 
IClaret fc Willemj (j2002|) revised the status of the apsidal 
motion test to stellar structure and evolution. They 
increased the observational sample by about 50% in 
comparison with previous works and took into account 
the effects of dynamic tides to determine the contribution 
of the tidal distortion to the predicted apsidal motion 
rate; they found a good agreement between observed and 
theoretical apsidal motion rates. The most recent internal 
structure co nstants for main-se q uence stars are those 
contained in IClaretl (|2004l . l2005l . l2006bl . l2007l hereafter 
C04, C05, C06b and C07, respectively). They were calcu- 
lated with new stellar models based on updated physics 
computed with different metallicities: (X, Z)=(0.70, 0.02) 
[C04]; (X,Z)=(0.754, 0.002) and (0.748,0.004) [C05]; 
(X,Z)=(0.730, 0.010) and (0.739,0.007) [C06b]; and 
(X,Z)=(0.64,0.04), (0.58,0.06) and (0.46,0.10) [C07]. 

2.1. Perturbations and the apsidal motion 

The longitude of periastron of a binary orbit, cu, defines the 
direction of the line of apsides in the orbital plane. It is 
constant, in the orbit of a system consisting of two gravi- 
tating bodies, only if all the three following conditions are 
valid: (i) the bodies can be regarded as point masses; (ii) 
they move in accordance with Newton law of gravitation 
(r~^); and (iii) the two bodies form a gravitationally iso- 
lated system. However, if any of these conditions fails, the 
size, form, and spatial position of the orbit will vary. The 
most readily detectable effect is a variation in the value of lo 
with time that is referred to as rotation (advance or reces- 
sion) of the line of apsides. For a more detailed discussion 
of this subject, see, fo r insta nce, the works of , Batten (19731 ) 
or IClaret fc GimenezI (|200lD . 

Several types of perturbations exist and can lead to ro- 
tation of apsides, such as mutual tidal distortion of the 
components, distortion of the components due to axial ro- 
tation, relativistic effects, presence of a third body, and re- 
cession due to a resisting circumbinary medium. The axial 
rotation and the mutual tidal forces of close binary sys- 
tems' components will deform each other and destroy their 
spherical symmetry, by means of the respective disturbing 
potentials. Besides the changes in the stellar structure, de- 
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scribed in Sect.[H these disturbing potentials produce an 
observed change in t o which is the sum of the effects caused 
by each component (jBattenl . [T97 3). The total rate of apsidal 
advance per orbital revolution, w, is 



2^ = c7 = ^'''-'^ 



(3) 



where P is the anomalistic orbital period, U is the apsidal 
motion period, and 



C2i 



fie) 



15M3 



Mi 



-.9(e) 



A . 



(4) 



where i=l,2 stands for the primary and the secondary stars, 
respectively; Mi and Ri are stellar mass and radius of com- 
ponent i; A is the semi-major axis; e is the orbital eccen- 
tricity; and the functions /(e) and g[e) are defined as 



/(e) = (1 - e2)-2 and .g(e) 



12e2 



2')/(e) 



2.5 



8 



(5) 



(rjj/wx) being the ratio between the actual angular rota- 
tional velocity of the stars and that corresponding to syn- 
chronization with the average orbital velocity. Equation ^ 
is a special case of Eq. ([2]), in which only the second or- 
der harmonics are taken into account. The first term in 
Eq. (HI) represents the contribution to the total apsidal mo- 
tion given by rotational distortions and the second term 
corresponds to the tidal contributions. 

With the exception of the k2i , all parameters in Eq. ^ 
can be determined from photometric and spectroscopic 
analysis. The empirical weighted average of the internal 
structure constants for individual systems can be given by 



1 



P 



1 



^2obs 



C21 + C22 U C21 + C22 2tt 



(6) 



Equations (jl]) and ([6]) show that fc2obs depend on 
our knowledge of the rotation velocities of the compo- 
nent stars. In most binaries with good absolute dimen- 
sions, the rotation velocities of the individual components 
are known through spectroscopic analysis. Since the av- 
erage orbital rotation, or Keplerian velocity, is a function 
of the orbital period, the ratio of rotational velocities in 
Eq. (|4l), namely Qi/ u jK, is well determined in these binaries 
(| Claret fc Gimenea . 119931 ). For systems without observa- 
tional determinations, the best approach is to assume that 
the components are synchronized with the orbital velocity 
at periastron, where the tidal for ces are at ma ximum. The 
rotation velocities are related by (jKopall . Il978[ ) 



(1 + e) .2 



(7) 



where wp is the angular velocity at periastron, e is the or- 
bital ec centricity, and ujk i s the K eplerian angular velocity 
(Eq.HD. IClaret fc Gimened ()1993f ) checked the vahdity of 
this approximation, achieving a good agreement between 
the observed and the predicted rotational velocities by as- 
suming synchronization at periastron (see their Fig. 6). 

The mean k2i values, obtained through Eq. ([6]), can 
be compared with those derived from theoretical models 
(Eq. fTU)) . However, the observed mean values of A:2obs should 
be first corrected from non-distortional effects, like rela- 
tivistic, third body, and interstellar medium contributions. 



2.2. Internal structure constants for spherically symmetric 
configurations 

Internal structure constants (kj) can be approximately 
computed based on the simple (and unrealistic) assump- 
tion that stars can be described by spherically symmet- 
ric models. These approximate kj have been computed for 
comparison with observed rates of aps idal motion. 

The Radau's differential equation (jKopall . Il959t ) is nu- 
merically integrated throughout all the structure with a 
4th-order Runge-Kutta method (jPress et al.l . ll99^ : 



dr 



■6^(r;,+l) + r7,(r,,-l)=j(j + l), 



(8) 



where 77j(0) = j — 2 (j=2,3,4), p{r) is the local density at 
a distance r from the center, and p{r) is the mean density 
within the inner sphere of radius r. The resulting value of 
the function r]j{R), which satisfies Radau's equation with 
R being the radius of the configuration, is used to obtain 



the individual values of fc. 



j + l-VjiR) 



(9) 



Because the observed motion of apsides is the sum of the 
motion produced by both stars, the quantities k2i (i=l,2 for 
the primary and the secondary star, respectively) cannot be 
obscrvationally determined separately. So, the theoretical 
counterpart of Eq. ^ is the weighted mean value of the 
individual theoretical apsidal motion constant. 



C2ifc: 



21,thoo 



+ C22fc22,theo 



«2,thco 



C2I + C22 



(10) 



where C21 and C22 are computed through Eq. ([4]) by using 
absolute dimensions, and fc2i,thoo and A;22,thco are the the- 
oretical apsidal motion constants for the primary and the 
secondary, respectively, obtained from the stellar models 
(Eq.[9]) for the corresponding mass and radius of each com- 
ponent. Before performing the apsidal motion comparison, 
one has to cheek if the models are able to reproduce basic 
stellar parameters, such as effective temperatures, and to 
predict a common age for the tw o components. A good dis - 
cussion of this subject is given bv lClaret fc Gimened (IT993I) . 
In Sect. ([7]) we show our predictions about apsidal motion 
for a chosen binary, comparing them with observational 
data. 



3. The Kippenhahn & Thomas formulation 

The KT70 method is a strategy for considering disturbing 
potentials in evolutionary stellar models more realistically, 
so that the distortion produced by a given disturbing po- 
tential is entirely included in the total potential function. 

To clarify how these disturbing effects were taken into 
account in the KT70 method, the equations are re-derived 
here. In this formulation, the spherically symmetric sur- 
faces, normally used in standard stellar models, are replaced 
by suitable non-spherical equipotential surfaces character- 
ized by the total potential tp, the mass enclosed by the 
corresponding equipotential surface whose surface area is 

and encloses a volume V^, and r^, the radius of the 
topologically equivalent sphere with the same volume l/J/,, 
enclosed by the equipotential surface. 
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For any quantity / varying over an cquipotcntial sur- 
face, we can define its mean value as 



(11) 



'<Aco„=t. 

wlrere 5*^, = ^ da, and da is the surface element. 

The local effective gravity is given by 5 = where 
dn is the (non-constant) separation between two successive 
equipotentials ip and ip + dTp, so that we have 



(9) = ^ I ^da, 



(»->^^ / (i) 



(12) 



(13) 



The volume between the surfaces ip and ip+dijj is given by 



dV[p= J dnda = dip 

lAconat. Ipcans 

from which we obtain 



[^)da^di;S^{g-'), (14) 



dip 



1 1 dA/^ 



(15) 



and the volume of the topologically equivalent sphere is 
given hy = ^r^. 

Eq. (jlSp can be combined with the general form of the 
hydrostatic equilibrium equation, 



dP 

dP _ GM^ 



to give 



dMj; 



/pi 



where /p is given by 



fp = 



47rr,4 1 



(16) 
(17) 

(18) 



With these corrections, the four stellar structure equations, 
with as the independent variable, become 



dP _ _GAUf 
^4 JP' 



dM^, 
dr^ 



dM^ 47rr^p' 

dL^, _ ^ ^dS 
dM^ ^ dt ' 

dT^ _ GM^.T _ ( ft 

dM^ 'iirr^P [ /p 

where fp is given by Eq. p8|) and 



47rr| 



1 



(19a) 
(19b) 
(19c) 
(19d) 

(20) 



In the case of isolated and non-rotating stars, /p=/t=l, 
and the original stellar structure equations are recovered. 
In order to obtain the internal structure of a distorted gas 
sphere, the set of Eqs. (fTO|) must be numerically integrated 
under suitable boundary conditions. 

This formulation was largely used in the literature 
mainly due to the easines s of its implementation in existing 



evolutionary codes (e.g. Endal fc Sofial. 19761: Lawl. 198C 



1999aD 



Pinsonneaull Il988t [Martin fc Claretl . Il996l : iMendes et al 



4. Tidal and/or rotational distortions on the 
equilibrium structure of stars 

We consider the effects of tidal forces, as well as the com- 
bined effects of tidal forces and rotation, on the stellar 
tructure and im plemented such effects in the ATON code. 
iMendesI ()l999bD had already introduced the effects of ro- 
tation alone in it. The tidal and rotational effects in stel- 
lar models were implemented accor ding to the treatmen t 
derived by KT70 and modified bv lEndal fc Sofial (|1976f ). 
who in turn used a more refined function to take into ac- 
count terms related to the distortion of the stars, namely 
the Clairaut-Legendre expansion for the grav itational po- 
tential of a self-gravitating body (jKopall . [19590 . The calcu- 
lations are done within the framework of static tides; in the 
case of dynamic tides, a more refined treatment is req uired 
(see iClaret fc Willed . [200I IWillems fc Claretl . [2001 . 

In our model, tides and axial rotation are the two phys- 
ical causes of the stellar configuration deviating from a 
spherical form. The centrifugal pseudo-potential terms aris- 
ing from the orbital motion around the center of mass of the 
system do contribute, also, for the departure from spheric- 
ity, but wc do not consider these terms in the present work, 
yet. These terms will be included in future works. We treat 
three different situations: 1) rotation acting alone, 2) tidal 
forces acting alone and 3) a combination of both effects dis- 
torting the star. The special case of the evolution of one of 
the components in a binary system is considered, assum- 
ing that (i) it rotates about an axis perpendicular to the 
orbital plane; (ii) from the orbital point of view, the sys- 
tem is considered to be at rest, i.e., the stars do not revolve 
around each other; and, particularly in this work, (iii) the 
ratio of the mutual separation to the radius of the evolving 
(distorted) star is assumed to be constant. The value of this 
ratio is an input parameter that depends on the system to 
be reproduced; in close binaries, the s eparation of the stars 
is often less than 10 times their radii ([Hilditchl . l200lf ). The 
cquipotcntial surfaces can be written as an expansion of the 
tesseral harmonics^ 



i,3 



(21) 



where ro is the mean radius of the corresponding cquipo- 
tential surface. The tesseral harmonics can be written as 

Yliro,9,cj^^K(ro)P;{e,<j>), (22) 

^ Tesseral harmonics are spherical harmonics of the form 
cos(m(^)P™(cos 9) and sin(m(j!>)P™(cos S), for m 7^ 1. These har- 
monics are so named because the curves on which they vanish 
are I — m parallels of latitude and 2m meridians, which divide 
the surface of a sphere into quadrangles whose angles are right 
angles. 
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where K{rQ) is the radial part of the tesseral harmonics 
and its angular part. P!, being functions associated with 
Legendre polynomials. 

4.1. Rotational distortion 

Rotation alone would render the star a rotational spheroid 
flattened at the poles. In what follows we are considering 
conservative rotation, which means that the centrifugal ac- 
celeration can be derived from a potential. 



(23) 



where s = r sin 9 is the perpendicular distance to the 
rotation axis of the star. A necessary and sufficient 
condition for the conse r vative case is that = 0(s) 
(jKippenhahn fc Weigertl . I1994D . which means rotation is 
constant on cylinders. Rigid body rotation is obviously a 

special case of conservativ e rotation.^ 

Following Eoial ([l959l) and later lEndal fc Sofial lfl97^ . 
the total potential is divided in three parts according to 
Eqs. (j24[) . where tps is the spherically symmetric part of 
the gravitational potential, ipr is the cylindrically symmet- 
ric potential due to rotation, and V'd™*^ is the cylindrically 
symmetric part of the gravitational potential due to the dis- 
tortion of the figure of the star caused by rotation. If the 
coordinates of the point P are the radius r and the polar an- 
gle 9 (measured from the rotational axis), the components 
of the potential at P can be written as 



1 



sin^ 6*, 



(24a) 
(24b) 



where we have abbreviated, as it is customary. 



m 



Yrot dro 



(29) 



In short, Yj-ot is a measure of the deviation from sphericity 
caused by rotation. The quantity 772 is of particular interest 
to our study, because the theoretical apsidal motion con- 
stant k2 can be derived from it (see Sect. [2]). The evaluation 
of 772 can be done by numerically integrating the Radau's 
equation: 



^— — r(?/2 + 1) + mim - 1) = 6. 



dro p{ro) 



(30) 



This equation is slightly different from Eq. ([8]). Here we use 
j = 2, and the spherical radius r was replaced by tq, the 
mean radius of the distorted configuration. 

By virtue of the definition of Krot = Y2 as given by Eq. 



, the expression for 7/'^'^°*'' simplifies to 



/ (rot) 



'3^ 



P2(C0S( 



V2 



M, 



dr' 



'4' 



V2 



(31) 



so that the total potential can be written as 



(rot) 



1 



sin^ 9 



ro y 



m 



ri2 



dr'o- 



(32) 



; (rot) 



00 

J— ^ n 



In Eqs. ([M)) . and Q. are respectively the mass and the 
angular velocity of the rotating star, while ro, by virtue of 
the cylindrically symmetric rotational potential ■0^, corre- 
sponds to the radius of the equipotential surface (as given 
by Eq. [21]) at the angle 9o defined such that P2(cos0o)=O, 
P2 being the second-order Legendre polynomial. Here, the 
shape of rotating configurations is described according to 
the expansion given in Eq. (|21[) , with given by 



Y;(ro) 



2j + l 



j + Vjiro)GM^ 



PU 



(25) 



By limiting ourselves to first-order theory, the disturb- 
ing potential associated to rotation (Eq. I24bp is such that 
it will invoke a single non-zero Y^ term in Eq. (|2ip , namely 
that one corresponding to i~0, .7 = 2 (jKopall . [19591 [l960l ) : 



By defining the radial part of the axisymmetric tesseral 
harmonic Y^ot as 



3GM^ 2 



(33) 



the radius of the equipotential surface, Eq. (|27p . can be 
rewritten as 



r{ro,9) = ro 1 - ^(ro)P2(cos( 



(34) 



To relate ro to r^, we evaluate the volume integral from 
r=0 to r(ro, 9) to obtain 



v., 



Anr^ 



_ 1^3 
5 35 



(35) 



For simplicity, the arguments of the term A(ro) were sup- 
pressed. From the previous equation, r^ is then obtained 
as 



C0,2 



(26) 



So, in the case of rotational distortions acting alone, the 
equipotential surface now becomes 



riro,9) = ro 



l + Y,r 



(27) 



in which we replaced Y2 by Yiot for better clarity, and 
5 



Y 



3GM^ 2 + 772(ro) 



P2(cos6l), 



r,p = ro 



1/3 



(36) 



Usually, r^ is known and so ro can be calculated through 
Eq. (j36p by means of an iterative procedure. 
Since the local effective gravity is given by 



(28) 9 



9-0 
dn 



dip 
dr 



1 dip 
r'd9 



1/2 



(37) 
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g can be found by differentiation of Eq. ([321) . The integral 
in Eq. p2p and their derivatives must be evaluated numer- 
ically. Once the values of (g) and {g~^) are known for a set 
of points on an equipotential surface, S^{g) and S-^{g~^) 
can be found, respectively, from Eqs. ([T^ and by nu- 
merically integrating over 9. 

4.2. Tidal distortion 

Tidal distortion acting alone would tend to elongate the 
star in the direction of the other component. Similarly to 
Eqs. fM]) . the total potential is divided in three parts ^s, "0*, 

and "^d*"^^ J where ips is again the spherically symmetric part 
of the gravitational, potential, ipt is the non-symmetric po- 
tential due to tidal forces, and V'd'"*^ is the non-symmetric 
part of the gravitational potential due to distortion of the 
star caused by the presence of the companion. If the coor- 
dinates of the point P are the radius r, the polar angle 6, 
and the azimuthal angle </>, the components of the potential 
at P can be written as: 



So, for tidal forces acting alone, the Yj are a measure of the 
deviation from sphericity due to those forces. 

As in Sect. l4.1i the apsidal motion constants can be de- 
rived from r]j by using our theoretical stellar models (see 
Sect. (2). The evaluation of the quantity above can be done 
by, again, numerically integrating the Radau's equation 



drn 



■(ry, +l)+r/,(r/j-l)=j(j + l), 



(44) 



for each j = 2, 3, 4. 

Substituting Yj as given by Eq. ((iT|) in Eq. (|38cp and 
carrying on the partial derivative in the latter, we get 



(tid) 



P,(A) 



(ri?)^+i 



J 



j + Vj 



'-Ildr'„. (45) 



The total potential is then 



(tid) 



r 

GM2 
R 



GM^ GM2 



f:(i)'p,(A) 

J=2 



(38a) 
(38b) 



R 

4 



Lj=2 



P,{X) 



R 



(ri?)^+i 
j-^ 



/2j+3 • , o 1 

(46) 



J+Vj 



By defining the radial parts of the non-symmetric 



^ AnG f d /j+3^iN , , /oo N By defining the radial part 
^ (2j -t- l)ri+^ J ^ dr' ^ ' °' ^ ' tesseral harmonics Y, in Eq. (|iT|) 

j = 2 r, 



In Eqs. ([38|) . M^, and A/2 are respectively the masses of the 
primary and the disturbing star; R is the mutual separation 
between the centers of mass of the two stars; X=cos(j)sin9; 
ro is the radius of the equipotential surface at the angles 
(^Oj^o), defined such that 



B(ro)P2(Ao) + C(ro)P3(Ao) + D(ro)P4(Ao) = 0, 



(39) 



M- 



M2 



M2_ /ro 
M^\R 

-(-) 



RJ 2 + ri2' 

4 7 



3 + %' 
5 9 



(47a) 
(47b) 
(47c) 



where Aq = cos 0o sin a-nd the terms B(ro), C(ro) and 
D(ro) arc defined in Eqs. (|47)) : and Pj are the jth-order 
Legendre polynomials. 

By virtue of the disturbing tidal potential given by Eq. 
(j38bp . and again restricting ourselves to first-order quan- 
tities only, it follows that the only non-zero Yj terms in 
Eq. (pij) consistent with the equilibrium theory of tides are 
those for wh ich i = 0, 7 = 2, 3. 4 (the reader is referred to 
the works of iKoDallfTOSl fT960l for fuU details): 



we can describe the equipotential surfaces, Eq. as 
r{ro,9,(f>) = 

ro [1 + B(ro)P2(A) + C(ro)P3(A) + D(ro)P4(A)] . (48) 

We evaluate the volume integral from r=0 to r{ro,0,(j)) 
given by Eq. (|48p and obtain 



AnrQ 



~5~ 



3^2 
7 



2P3 
"35" 



Co 



G- 



M2 



PJ+i 

from which we have 
M2 2j + 1 



(j = 2,3,4), 



rn \ 



j + Viiro) 
where, as before, 

"^'^Y.dro' 

the equipotential surfaces are then described by 

4 



(40) 



(j= 2,3,4), (41) 



35 

6C^D 18L»3 



3 



35 



231 



77 



1001 



(49) 



where the arguments of the terms B{ro), C{ro) and D{ro) 
were omitted. From Eq. (j49p . rg and r^ are related by 



ri, = ro 



3P2 



7 



2P3 
~35~ 



6B^D 4:BC^ 



35 



35 



riro, 



ro 



1 



(42) 



(43) 



+- 



231 



77 



1001 



3 



(50) 



Similarly to the case of pure rotation, tq can be calculated 
from the value of through Eq. (j50p , by means of an iter- 
ative procedure. 
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The local effective gravity is now given by 



Il989l ). the external equipotential surface of the distorted 
star can be described by 



dn 



dt/j \ ^ /I dt/j \ ^ / 1 9-0 



dr 



86 



r sm ( 



(51) 



(54) 



and it can be obtained by differentiation of Eq. . As 
before, integrals and derivatives must be evaluated numer- 
ically. 

Finally, by integrating over 9 and </>, we find S.,p{g) and 
Sjp{g~^) from Eqs. (fT2|) and p3)l . respectively, once (g) and 
{g~^) are known for a set of points on an equipotential 
surface. 

4.3. Combined effects of rotation and tides 



According to iKopall (|l960l Il974l ) , the combined effects of 
rotation and tides give rise to a total disturbing potential r(ro,6', (/>) = 
that can be expressed as the sum of the distortion terms 
corresponding to pure rotation or tides and of a number 
of interaction terms. The latter, however, are in general 
second-order terms; so, in the present work, we neglect 
those interaction terms since we are dealing with first-order 
theory only. 

The total potential for the combined case is then com- 
posed by four parts, which represent the joint contributions 
from first-order rotational and tidal effects as respectively 
given by Eqs. ([M]) and ([351): ips, the spherically symmet- 
ric part of the gravitational potential; ipr, the cylindrically 
symmetric potential due to rotation; ipt, the non-symmetric 

potential due to tidal forces; and -0^ = V'd™*'' + the 
non-symmetric part of the gravitational potential due the 
to distortion of the star considering both effects. Again, 
we do not take into account the pseudo-potential centrifu- 
gal terms due to the orbital motion in this approximation. 
So, the total potential for the combined case (to first-order 
approximation) at P{r, 9, 0) is 



where y,ot is given by Eq. ((28)) and the Yj are given by 
Eq. (gl]). As in Sects. (|i?T|) and (g^l), the evaluation of rjj 
can be done by numerically integrating Radau's equation 
fEa. l44p . In this way, by using an adequate iterative proce- 
dure, we can obtain a value for 772 that reflects the combined 
effects of rotation and tides, and 773 and 774 values that ob- 
viously relate only to tidal effects. 

If we define the radial parts of the tesseral harmonics, 
A(ro), B{ro), C(ro) and D{ro), according to Eqs. ([55)) and 
([TT)) . Eq. ([5l)) can be rewritten as 



ro[ l-A(ro)P2(cos0) + B(ro)P2(A) 

+C(ro)P3(A)+i^(ro)P4(A) ]. (55) 

By integrating Eq. (|55p from r=0 to r(ro, 9, 0), we obtain 

SAB 



3^2 

lOAD^ 

231 
6C^D 



140 ^ 
QABD 

20BD^ 
231 



3^2 3^2 2^3 
~5~ ^ ~ 35" 
'SAB^ 2AC^ 2B^ 
^ ^ ~35~ 



35 
6B^D 

' 35 
IBL'^ 



35 
4BC2 



35 



1001 



3 



(56) 



For simplicity, A(ro), -B(ro), C(ro), D{rQ) appear with no 
arguments. From Eq. ()56p . tq and are related as 



i' 



Ips + i'r +1pt+1pd 



—ir sm 9 H 

2 R 



1 + E 



i=2 



R 



P,{X) 



7 

4 

-47rGA/2 ^ 



3^2 


iAB 


3B^ 


3C2 2^3 




+— - 


' 5 + 


— + 


~ 35" 




iA^B 


\ 


3AB'^ 


2ylC2 




2^3 


35 ' 


140 


35 


+ 35 




~35~ 


iA^B 


QA^D 
\ 


3AB^ 


2AC^ 


+ 


253 


35 ' 


140 


35 


+ 35 








20SD2 




Z?2 1 


1 
3 




77 


^ 231 


^ 1001 


+ 3 







0=2 



i + Vj 



dr'o. (52) 



In Eq. (|52p . M^, A/2, ^, R, and A retain their meanings as 
previously defined in Sects. ()4.ip and ()4.2p . while 7'o is now 
the radius of the equipotential surface at the angles (^oi 0o) 
defined such that 



-A(ro)P2(cos0o) + B(ro)P2(Ao) + C(7'o)P3(Ao) 
+D(ro)P4(Ao) = 0, 



(53) 



where A(ro) is given by Eq. ([55)) . and 6(7-9), C(ro), and 
D(ro) are defined in Eqs. ([47)) . 

Since, as previously mentioned, the distortions arising 
from both rotation and tidal effects are simply addit i ve (to 
first- and even second-order accuracy; lKopall[T960l . [197I 



(57) 



As before, tq can be calculated through Eq. ([57)) by means 
of an iterative procedure. 

The local effective gravity is given by Eq. ([5T)) . g can 
be numerically found by differentiation of Eq. (|52p . With 
the values of (g) and {g~^) known for a set of points on 
an equipotential surface, S^{g) and S*^, (5"^) can be found, 
respectively, from Eqs. ()12p and ([T5)) by numerically inte- 
grating over 9 and 0. 

4.4. Rotational inertia 

Rotational inertia (or moment of inertia) is an important 
tool for studying tidal evolution theories, since it is re- 
quired to p redict the c ircularization and synchronization 
ti me scal e s (IZahnl . Il977[ ). 

iMotzl ([l952f ) computed values of rotational inertia by 
using the existing (and unrealistic) models at that time. 
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iRucinskil ()1988f ) obtained values of the radii of gyration 
for polytropic models of low-mass stars at the ZAMS. The 
radius of gyration of a body is the distance between a given 
axis of this body^ and its center of gyration'^, and is defined 

by 



(58) 



where / is the rotational inertia of the star, and M 
and R are the s tellar mass and radius, respectively. 
IClaret fc GimeneS (jl989bl ) presented radii of gyration cal- 
culations for more massive stars during the hydrogen burn- 
ing phases. They used standard models, in which the 
stars are described by spherically symmetric configurations. 
More recent computations of radii of gyration were pro- 
vided by C04, COS, C06b, and COT. 

As a consequence of the distortions introduced by rota- 
tion and tides, the rotational inertia of the star is cha nged 
from that one corresponding to spherical symmetry. iLawl 
(jl980t ) derived the rotational inertia of a rotationally dis- 
torted mass shell, to first-order accuracy, as 



AI = -dm. 



- , ^0 

^^^^ ^ ^ 



where 



Pa^lcose*) sin^ 9(19. 



(59) 



(60) 



In what follows, we present a new expression for the 
rotational inertia taking into account both rotational and 
tidal distortions, again considering only first-order effects. 
We start by considering a very thin mass shell, for which 
the rotational inertia is given by 



AI = / RMm. 



(61) 



27r 7T r2 



AI ^ p 



r"* sin^ 9drd9d(l) 



ri 
27r 7r 



{rl - rf) sin^ 



(62) 



where ri and r2 are respectively the inner and outer radius 
of the mass shell. For a tidally and rotationally distorted 



In a rotating body, the rotation axis is considered. If no 
axis is specified, the centroidal axis, which is the line joining 
the centroid of each cross section along the length of an axial 
member such as truss diagonal, is assumed. 

^ The center of gyration of a body is defined as that point 
at which the whole mass might be concentrated (theoretically) 
without altering the body's rotational inertia. In other words, 
this is the center about which the body can rotate without mov- 
ing linearly or vibrating. 



mass shell, ri and r2 correspond to the radius of the equipo- 
tential surfaces of the distorted configuration as given by 
Eq. inil), so that 



(63) 



ro,[ I- A{ro^)P2{cose)+B{ro,)P2{\) 
+ C{ro,)P-i(\) + D{ro,)Pi{\) 1, 



where i=l,2. 

Since the terms B{rQi), C(roi) and D{rQi), in Eq. ([55]) . 
are proportional to {roi/Ry with j = 2,3,4 respectively, 
and that, furthermore, they would be raised to the fifth 
power when substituted in Eq. ([5^ . it follows that B{roi) 
would dominate over C(roi) and D{rQi) in the latter equa- 
tion, as rQi/R< 1. This, in short, means that terms of order 
j higher than 2 do not contribute significantly to the total 
departure from the spher ical symmetry, as already noted by 
IClaret fc Willemd (I2002D . So, as a further approximation, 
we drop the terms corresponding to j = 3,4 in Eq. (j63p 
for the sole purpose of computing the rotational inertia of 
a mass shell distorted by tides and rotation. This approxi- 
mation allows us to reduce from 252 to 42 the terms that 
result from the full development of Eq. ([5^ when ri and 
r2 are substituted for those given by Eq. ([63]) . In this way, 
AI becomes, after some algebra. 



AI 



Jo 



J2i-ly4^[l-A{r0^)P2{cO8 9) + 
i=l 

+B{rM)P2{X)^^\sm^0d0d(l). (64) 



The fifth power term within brackets in Eq. (|64p can be 
expanded by using the multinomial theorem 

{Xi + X2 + ... 4- !>)" ^ 

n 

, oi, a2, Op 



E 

2+ ■- + 

yielding 



0<ai,a2,.-.,ap<n 



"1 -^2 



(65) 



By using spherical coordinates, so that R = j'sin^ and 
dm = pr^ sm9drd6d(j), and assuming that p is constant in /\/; 
such a thin mass shell (which will be later justified as we 
get to the limiting differential case), the previous expression 
becomes 



0<ai,a2,a3<5 
01+02 + 113 = 5 



^(-l)*+Vo^^"^(roOB"n 



(66) 



where 



Qa=\ , / / P^'{cos9)P^^X)sm'' 9d9dcl). (67) 

\yai\a2\a3\ J Jq Jq 

For the limiting case of a thin mass shell, we have 

roi = ro, 

ro2 = roi+ dr 01 = tq + dro, 

A{roi) = A, 

B{roi) = B, 

A{ro2) = A{roi) + dA{rQi) = A + dA and 

B{ro2) = B{roi) + dB{roi)=B + dB, 

so that the terms within brackets in Eq. (|66p can be rewrit- 
ten as 



-rl2A^- {ro2)B-'' (ro2) + r^.A^' (roi )B"-^ (roi ) 



-(ro + dro^iA + dA^' {B + dBf 



.^5^a2^a3 



• (68) 
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For a thin mass shell one has j'q ^ dro and consequently 
A:$>dA and B > dB, so the right side of Eq. dHH]) can be 
further reduced to 



(rot) 



(tid) 



where 772'°*'' and 772*"^'' stand for the 772 values that would 
result from pure rotation or pure tidal effects, respectively. 

By using this latter result and remembering that dm^ = 
Airpr^dr^, Eq. ([66)1 can be brought to the form 



AI 



zOn rl dr. 



ro dro 



E 



(rot) 
02% 



(tid) 

asm 



(69) 



where Qa is given by Eq. ((67)) . 

In the special situation in which only tidal forces are 
present, the rotational inertia of a given mass shell is 



AT "^.1 2froYdro 
AI=-dm,r,(-) — 



1 + ^EQp^'(^^^2 + 5) 

Jtid) 



, (70) 



where 772 corresponds unambiguously to and Qp is 

given by 



p!(5-p)!yo 



27T n-K 



PP{X) sin^ 9d9d(t). 



(71) 



stars and finally (iv) rotating binary models (distorted si- 
multaneously by rotation and tidal forces), useful to study 
rotating stars in binary systems. In cases in which rotation 
is present, we assume rigid body rotation. The relation be- 
tween initial angular momentum (Jin) and stellar mass that 
we used was obtained from the respect ive mass-radius and 
mass- moment of inertia relations from iKawalerl ()1987f ): 

Jkaw - 1.566 X lO"^" i^—j cgs. (72) 

For computing the binary models (both rotating and non- 
rotating) we assumed a separation of 7 times the radius of 
the star whose evolution is followed (from this point on re- 
ferred to as the prima ry), which is in the range typical for 
close binary systems ()Hilditchl . |2001[ ) ; the disturbing star 
is considered to be a point mass of the same mass as its 
primary. As previously stated, the pseudo-potential cen- 
trifugal terms arising from the orbital motion around the 
system's center of mass are not included in the present ap- 
proximation for the tidal effects due to a companion. We 
followed the evolution of internal structure constants and 
moment of inertia during the pre-MS phase and tabulated 
them together with the corresponding evolutionary tracks. 
Table[l] presents the 1 Mq standard model as an example 
of such tables. Column 1 gives the logarithm of stellar age 
(in years); column 2, the logarithm of stellar luminosity (in 
solar units); column 3, the logarithm of effective tempera- 
ture (in K); and column 4, the logarithm of effective gravity 
(in cgs). Columns 5, 6 and 7 give the logarithm of internal 
structure constants, and column 8 gives the gyration radius. 



5. Models 

Our new version of the ATON evolutionary code has not 
only many updated and modern features, reg arding the 
physics of stellar interior (| Ventura et al.l . Il998[ ). but, also, 
is able to reproduce stars with spherically symmetric con- 
figurations, as well as tidally and rotationally distorted 
stars (see Sect.]?]). The grids cover a mass range from 
0.09 to 3.8 Mq and were computed from early stages 
of pre-MS phase up to the m ain sequence. Th e radia- 
tive opacitie s are taken from Iglesias fc RogersI (|1993() . 
extended by [Alexander fc FergusonI ( "994 ) tables in the 
low-temperature reg ime. The OPAL equation of state 
([Rogers et al.l . Il996[) is used in the range 3.7<logT<8.7, 
while in the low -T high density regim e we use the equa- 
tion of state of iMihalas et al.l (|1988D . The nuclear net- 
work includes 14 elements an d 22 reactions; the r elevan t 
cross-sections are taken from ICaughlan fc Fowled ()1988f ). 
We adopted the solar metalhcity (Z=0.017 5 y=0.27). The 
classi cal Mixing Length Theory (MLT, iBohm-Vitensd . 
[19581) was used to treat the convective transport of energy, 
though the ATON code can also use the EST (Fu ll Spectrum 
of Turbulence) treatment l|Canuto et al.l . [T996[ ). The mix- 
ing length parameter has been set to a=A/ Hp=1.5, value 
that, according to our calibration, best reproduces the solar 
radius at the solar age with b oundary conditions o btained 
from gray atmosphere models (|Landin et al.l . [20061 ). 

We present four sets of evolutionary models, namely: 
(i) standard, spherical models (with no distorting effect) 
corresponding to single non-rotating stars, (ii) binary mod- 
els (distorted only by tidal forces) in which we treat non- 
rotating stars in binary systems, (iii) rotating models (dis- 
torted only by rotation) for representing single rotating 



Table 1: Pre-MS evolutionary tracks (including log kj and 
(3) for IM0 star generated with our standard models". 



Jg Age 

(yr=) 



log - 



(K) 



off 



log g 



log fc2 log fca 



log k4 



2.2954 
3.6550 
4.1464 
4.5435 
5.1029 
5.3018 
5.5216 
5.7953 
6.0918 
6.3961 
6.7156 
7.0039 
7.2246 
7.3696 
7.4906 
8.7585 
9.3538 
9.5876 
9.6979 
9.8058 



1.8909 
1.6727 
1.4526 
1.2322 
1.0251 
0.8982 
0.6905 
0.4699 
0.2493 
0.0286 
-0.1866 
-0.3157 
-0.2773 
-0.1155 
-0.0389 
-0.1301 
-0.0815 
-0.0307 
0.0084 
0.0619 



3.5902 
3.6057 
3.6185 
3.6293 
3.6376 
3.6418 
3.6472 
3.6510 
3.6526 
3.6516 
3.6497 
3.6560 
3.6840 
3.7255 
3.7573 
3.7516 
3.7553 
3.7587 
3.7608 
3.7627 



1.8604 
2.1405 
2.4119 
2.6755 
2.9159 
3.0597 
3.2889 
3.5245 
3.7517 
3.9682 
4.1757 
4.3301 
4.4038 
4.4078 
4.4585 
4.5270 
4.4933 
4.4558 
4.4252 
4.3793 



-0.9075 
-0.8654 
-0.8325 
-0.8060 
-0.7868 
-0.7776 
-0.7664 
-0.7587 
-0.7535 
-0.7769 
-0.7904 
-0.8492 



0944 
5590 
8961 
7105 
8011 
9017 
9810 



-2.0934 



-1.3965 
-1.3526 
-1.3076 
-1.2707 
-1.2434 
-1.2301 
-1.2139 
-1.2028 
-1.1954 
-1.2180 
-1.2420 
-1.3026 
-1.5458 
-2.0471 
-2.4599 
-2.2631 
-2.3570 
-2.4620 
-2.5432 
-2.6579 



-1.7387 
-1.6927 
-1.6379 
-1.5920 
-1.5574 
-1.5405 
-1.5197 
-1.5053 
-1.4958 
-1.5174 
-1.5475 
-1.6126 
-1.8532 
-2.3655 
-2.8086 
-2.6077 
-2.7027 
-2.8096 
-2.8916 
-3.0078 



0.4177 
0.4253 
0.4319 
0.4374 
0.4415 
0.4435 
0.4459 
0.4476 
0.4488 
0.4495 
0.4458 
0.4260 
0.3811 
0.3202 
0.2812 
0.2974 
0.2872 
0.2767 
0.2689 
0.2588 



'^Thc complete version of the table, including 88 tracks for the four 
sets of models and all the masses of Tablef2l will be available only 
in electronic form. The tracks corresponding to rotating models and 
rotating binary models will contain an additional column regarding 
the rotational period. 



Table [2] presents the internal structure constants and 
the gyration radius at the ZAMS, for both non- rotating 
and rotating models. Its successive columns give the stellar 
mass (in Mq)] the logarithm of the stellar luminosity (in 
solar units); the logarithm of the effective temperature (in 
K); the logarithm of the surface gravity (in cgs units); the 
logarithm of the internal structure constants and 
A:4; the radius of gyration (3 (cf. Eq. [58[) : and, finally, the 
rotational periods (in days), when rotation is present. 



Landin et al.: Tidal and/or rotational distortions 11 



Table 2: Internal strueture eonstants and gyration radii for ZAMS rotating models. See the text for details. 



M 
(Mq) 


lo 




log ^cff 

(K) 




log g 
(cgs) 


log k2 


log fca 


log k4 




p 




p 

(d) 


M 
(Mg) 


lo 




log J ott 




logs 


log 


1 


ogfcs 


log ki 




p 




p 

(d) 






standard (s 


pherically symmetric) models 


















tidally 


(non-rotating) distorted models 










0.09 


-3 


3197 


3.4367 


5 


4113 


— 0.8134 


_r 


2455 


— 1.5272 





4600 









09 


-3 


3260 


3.4355 


5 


4146 


— 


8173 


_ 1 


2500 


— \ 


5318 





4601 






0.10 


-3 


0098 


3.4839 


5 


3360 


— n 8359 


_r 


2743 


— 1.5613 





4554 









10 


-3 


0131 


3.4836 


5 


3398 


— 


8391 




2779 


— \ 


5651 





4556 






0.20 


-2 


2529 


S 5255 


5 


0462 


— 0.8943 


_r 


3520 


— 1.6567 





4453 









20 


— 2 


2540 


3.5255 




0490 


— 


9000 


_ \ 


3592 


_ \ 


6651 





4454 






30 


— 1 


9334 


3 5444 


4 


9788 


— 8720 


_4 


3253 


— 1 6269 





4489 









30 


— 1 


9344 


3 5445 


4 


9815 


— 


8770 


_ \ 


3317 


_\ 


6345 





4490 






0.40 


-1 


6724 


3^5626 


4 


9153 


-0^9176 


_i 


3812 


— 1^6922 





4402 









40 


-1 


6716 


3.5628 


4 


9168 


— 


9213 


_1 


3870 


_1 


6995 





4397 






0.50 


— 1 


3820 


3^5885 


4 


8254 


-1^0460 


-1 


5023 


-1^8099 





4119 









50 


— 1 


3815 


3.5888 


4 


8275 


-1 


0501 


-1 


5082 


-1 


8173 





4114 






0.60 


— 1 


0632 


3.6252 


4 


7325 


-1.2494 


-1 


7065 


-2.0127 





3808 









60 


— 1 


0598 


3.6240 


4 


7260 


-1 


2309 


-1 


6897 


-1 


9975 





3821 






0.70 


— 


7584 


3.6639 


4 


6496 


-1.4603 


-1 


9300 


-2.2404 





3534 









70 


— 


7558 


3.6640 


4 


6487 


-1 


4621 


-1 


9334 


-2 


2449 





3532 






0.80 


-0 


4808 


3.7026 


4 


5849 


-1.6441 


-2 


1288 


-2.4442 





3314 









80 


-0 


4778 


3.7024 


4 


5825 


-1 


6417 


-2 


1282 


-2 


4453 





3312 






0.90 


-0 


3454 


3.7280 


4 


6022 


-1.7067 


-2 


2259 


-2.5554 





3216 




— 





90 


-0 


2348 


3.7312 


4 


5059 


-1 


8397 


-2 


3484 


-2 


6735 





3071 




— 


1.00 


-0 


1253 


3.7540 


4 


5319 


-1.9037 


-2 


4571 


-2.8011 





2957 






1 


00 


-0 


1228 


3.7546 


4 


5331 


-1 


9185 


-2 


4769 


-2 


8234 





2950 






1.20 





3184 


3.7935 


4 


3251 


-2.4573 


-3 


1435 


-3.5662 





2361 






1 


20 





3319 


3.7941 


4 


3156 


-2 


4968 


-3 


1894 


-3 


6146 





2334 






1.40 





6830 


3.8214 


4 


1394 


-2.8859 


-3 


7014 


-4.2406 





1962 






1 


40 





6901 


3.8219 


4 


1355 


-2 


9077 


-3 


7278 


-4 


2709 





1950 







1.60 





9558 


3.8590 


4 


0749 


-2.9992 


-3 


8561 


-4.4404 





1861 






1 


60 





9630 


3.8595 


4 


0709 


-3 


0239 


-3 


8827 


-4 


4670 





1853 







1.80 


1 


1839 


3.8998 


4 


0611 


-2.9891 


-3 


8436 


-4.4275 





1869 






1 


80 


1 


1883 


3.9011 


4 


0629 


-3 


0035 


-3 


8618 


-4 


4467 





1866 







2.00 


1 


3797 


3.9357 


4 


0547 


-2.9667 


-3 


8203 


-4.4044 





1887 






2 


00 


1 


3851 


3.9366 


4 


0541 


-2 


9796 


-3 


8339 


-4 


4182 





1882 






2.30 


1 


6345 


3.9811 


4 


0419 


-2.9337 


-3 


7823 


-4.3644 





1913 






2 


30 


1 


6376 


3.9826 


4 


0464 


-2 


9370 


-3 


7882 


-4 


3727 





1912 







2.50 


1 


7839 


4.0073 


4 


0339 


-2.9098 


-3 


7554 


-4.3376 





1930 






2 


50 


1 


7867 


4.0089 


4 


0385 


-2 


9233 


-3 


7718 


-4 


3542 





1930 







2.80 


1 


9782 


4.0438 


4 


0349 


-2.8593 


-3 


7029 


-4.2840 





1967 






2 


80 


1 


9834 


4.0447 


4 


0343 


-2 


8732 


-3 


7160 


-4 


2976 





1962 






3.00 


2 


0976 


4.0649 


4 


0296 


-2.8407 


-3 


6809 


-4.2612 





1985 






3 


00 


2 


1002 


4.0664 


4 


0344 


-2 


8480 


-3 


6905 


-4 


2709 





1985 






3.30 


2 


2587 


4.0938 


4 


0257 


-2.8049 


-3 


6406 


-4.2186 





2014 






3 


30 


2 


2589 


4.0959 


4 


0350 


-2 


7987 


-3 


6360 


-4 


2145 





2018 







3.50 


2 


3542 


4.1121 


4 


0287 


-2.7741 


-3 


6065 


-4.1826 





2037 






3 


50 


2 


3593 


4.1128 


4 


0275 


-2 


7974 


-3 


6341 


-4 


2126 





2032 







3.80 


2 


4869 


4.1368 


4 


0307 


-2.7394 


-3 


5731 


-4.1496 





2069 






3 


80 


2 


4915 


4.1378 


4 


0313 


-2 


7537 


-3 


5871 


-4 


1636 





2065 






rotationally distorted (single star) models 










rotationally an 


d tidally distorted models 










0.09 


-3 


5403 


3.3830 


5 


4499 


-0.8321 


-1 


2660 


-1.5486 





4657 





0323 





09 


-3 


5459 


3.3818 


5 


4520 


-0 


8317 


-1 


2657 


-1 


5481 





4655 





0322 


0.10 


-3 


2999 


3.4264 


5 


4644 


-0.8653 


-1 


3053 


-1.5923 





4666 





0237 





10 


-3 


3060 


3.4251 


5 


4666 


-0 


8649 


-1 


3047 


-1 


5916 





4664 





0236 


0.20 


-2 


3178 


3.5083 


5 


1457 


-0.9989 


-1 


4770 


-1.7971 





4450 





0430 





20 


— 2 


3170 


3.5086 


5 


1456 


-0 


9965 


-1 


4741 


-1 


7938 





4449 





0429 


0.30 


-1 


9701 


3.5350 


5 


0252 


-0.9202 


-1 


3823 


-1.6909 





4520 





0777 





30 


— 1 


9707 


3.5349 


5 


0269 


-0 


9194 


-1 


3815 


-1 


6901 





4518 





0775 


0.40 


-1 


7068 


3.5551 


4 


9499 


-0.9340 


-1 


4047 


-1.7203 





4465 





1159 





40 


— 1 


7041 


3.5553 


4 


9491 


-0 


9358 


-1 


4072 


-1 


7232 





4457 





1160 


0.50 


-1 


4194 


3.5808 


4 


8578 


-1.0457 


-1 


5085 


-1.8213 





4182 





1561 





50 


— 1 


4171 


3.5810 


4 


8577 


-1 


0438 


-1 


5063 


-1 


8189 





4176 





1562 


0.60 


-1 


1005 


3.6164 


4 


7590 


-1.2446 


-1 


7069 


-2.0180 





3852 





1995 





60 


— 1 


0981 


3.6165 


4 


7583 


-1 


2424 


-1 


7044 


-2 


0152 





3849 





2001 


0.70 


-0 


7913 


3.6551 


4 


6710 


-1.4587 


-1 


9288 


-2.2416 





3566 





2445 


n 
u 


/ u 


-0 


7897 


3.6550 


4 


6706 


-1 


4540 


-1 


9246 


-2 


2378 





3565 





2454 


0.80 


-0 


5102 


3.6945 


4 


6054 


-1.6515 


-2 


1405 


-2.4603 





3334 





2845 





80 


-0 


5075 


3.6935 


4 


6000 


-1 


6403 


-2 


1236 


-2 


4405 





3339 





2897 


0.90 


-0 


2657 


3.7247 


4 


5344 


-1.8640 


-2 


3763 


-2.7039 





3084 





3243 





90 


-0 


2638 


3.7248 


4 


5342 


-1 


8639 


-2 


3754 


-2 


7025 





3080 





3247 


1.00 


-0 


1458 


3.7465 


4 


5486 


-1.9242 


-2 


4777 


-2.8225 





2963 





3233 


1 


00 


-0 


1444 


3.7467 


4 


5495 


-1 


9231 


-2 


4757 


-2 


8197 





2959 





3228 


1.20 





2881 


3.7838 


4 


3569 


-2.4765 


-3 


1529 


-3.5633 





2362 





3968 


1 


20 





3355 


3.7672 


4 


2349 


-2 


4134 


-3 


0060 


-3 


3663 





2432 





5464 


1.40 





6679 


3.8070 


4 


1585 


-2.9872 


-3 


8173 


-4.3553 





1891 





4931 


1 


40 





6753 


3.8065 


4 


1514 


-3 


0084 


-3 


8409 


-4 


3800 





1878 





4974 


1.60 





9437 


3.8406 


4 


0753 


-3.1595 


-4 


0570 


-4.6681 





1749 





5857 


1 


60 





9492 


3.8400 


4 


0696 


-3 


1646 


-4 


0617 


-4 


6729 





1740 





5902 


1.80 


1 


1741 


3.8852 


4 


0602 


-3.1158 


-4 


0057 


-4.6137 





1778 





6840 


1 


80 


1 


1772 


3.8852 


4 


0587 


-3 


1232 


-4 


0129 


-4 


6212 





1774 





6854 


2.00 


1 


3757 


3.9229 


4 


0452 


-3.0779 


-3 


9557 


-4.5578 





1807 





7948 


2 


00 


1 


3793 


3.9225 


4 


0415 


-3 


0858 


-3 


9640 


-4 


5662 





1801 





7992 


2.30 


1 


6297 


3.9723 


4 


0388 


-3.0079 


-3 


8755 


-4.4723 





1858 





9586 


2 


30 


1 


6333 


3.9719 


4 


0353 


-3 


0205 


-3 


8886 


-4 


4851 





1852 





9638 


2.50 


1 


7790 


4.0003 


4 


0331 


-2.9768 


-3 


8396 


-4.4328 





1886 


1 


0773 


2 


50 


1 


7806 


4.0007 


4 


0342 


-2 


9814 


-3 


8443 


-4 


4376 





1884 


1 


0762 


2.80 


1 


9768 


4.0378 


4 


0290 


-2.9217 


-3 


7756 


-4.3644 





1927 


1 


2590 


2 


80 


1 


9787 


4.0380 


4 


0296 


-2 


9249 


-3 


7803 


-4 


3700 





1925 


1 


2588 


3.00 


2 


0944 


4.0602 


4 


0286 


-2.8859 


-3 


7392 


-4.3279 





1954 


1 


3821 


3 


00 


2 


0957 


4.0607 


4 


0309 


-2 


8819 


-3 


7328 


-4 


3199 





1954 


1 


3785 


3.30 


2 


2542 


4.0906 


4 


0287 


-2.8342 


-3 


6800 


-4.2646 





1993 


1 


5722 


3 


30 


2 


2577 


4.0904 


4 


0260 


-2 


8433 


-3 


6885 


-4 


2723 





1988 


1 


5784 


3.50 


2 


3530 


4.1085 


4 


0259 


-2.8139 


-3 


6572 


-4.2405 





2014 


1 


7091 


3 


50 


2 


3527 


4.1094 


4 


0311 


-2 


8084 


-3 


6511 


-4 


2333 





2016 


1 


6983 


3.80 


2 


4866 


4.1338 


4 


0271 


-2.7748 


-3 


6134 


-4.1936 





2048 


1 


9082 


3 


80 


2 


4880 


4.1343 


4 


0291 


-2 


7711 


-3 


6083 


-4 


1878 





2047 


1 


9038 



6. Results and discussion 

In this section, we discuss the results yielded by the models 
presented in Sect. ([5]) and the differences obtained with each 
set of models. Then we proceed with a comparison with 
other works in the literature. 

6.1. Effects on the internal structure constants 

As previously mentioned, iMohan et al.l ljl990[ ) found that 
the effects of rotational distortions on main-sequence stellar 
models are greater than those of tidal distortions, though 
their results were obtained only for binary models with a 
mass ratio of 0.1. We investigated the ZAMS models in or- 
der to verify such differences, by calculating how different 
the non-standard values of the second-order internal struc- 
ture constants are as compared with the standard ones for 
each stellar mass at the ZAMS. The binary models pro- 
duced values of logfc2 lower than the standard models, on 
average, by a difference of about 0.0170, with a maximum 
difference of 0.1330 occurring for the 0.9 M© model. For 
the rotating and the rotating binary models this average 
difference is 0.0570 and 0.0594, respectively; the maximum 
difference happens at the mass of 1.6 Mq for both sets of 
models, being 0.1603 for rotating models and 0.1654 for 
rotating binary ones. 

To better assess the relative importance of the rota- 
tional and tidal effects, we also ran rotating binary models 
of 0.3 and 1.0 solar masses with (a) separations of 3, 14, 35 
and 70 times the radius of the primary, (b) masses of the 



disturbing star corresponding to the double and half val- 
ues of the primary mass, and (c) initial angular momentum 
equal to 0, 1 and 3 Jkaw, where Jkaw is given in Eq. ([7^ . 
Table[3] presents, from left to right, the primary mass, the 
initial angular momentum, the secondary mass, the separa- 
tion between the components and the logarithm of ^2 at the 
ZAMS for these additional models. From that table we can 
see that, with the only exception of the case of 1 Mq with 
orbital separation of 3 times the primary's radius and Jjn 
equal to Jkaw, logA:2 is mainly affected by the inclusion of 
rotational effects. For 0.3 M©, considering wider orbits and 
assuming different mass binary components do not alter 
significantly the log ^2 values. For the 1 Mq model, log ^2 
decreases slightly when the distance between the stars in- 
creases, except for separations of 3 and 70 stellar radius, 
and when we vary the mass of the companion stars. The 
maximum difference occurs for 1 Mq models, when we con- 
sider a companion star of 2 Mq . 

The particular cases of 1 M© with orbital separation of 
3 times the primary's radius and Jin = Jkaw or Jin — 3 Jkaw 
(where Jkaw denotes the initial angular momentum given by 
Eq.[75]) demonstrate however that, depending on the star's 
rotation rate, a shorter orbital separation can make the 
tidal effects to overcome the rotational ones regarding the 
internal structure constants. This can be seen in Fig. ([T|), 
which also shows that the differences become more apparent 
for ages starting at the ZAMS. 

For the gyration radii, we found, for a given mass at 
the ZAMS, values slightly higher for the rotating binary 
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log Age (Myr) 



Fig. 1: logfc2 values against age for an IMq models with orbital separation R and initial angular momentum Jjn as 
indicated, where r^, denotes the primary's radius and Jkaw the initial angular momentum given by Kawaler's l|l987f) 
relations. Long-dashed lines represent values for the standard (non-rotating, single star) model; solid lines, the rotating- 
only model; dotted lines, the binary model without rotation; and the short-dashed lines, the rotating binary model. 



Table 3: Values of logfc2 obtained for additional models'". 



Mprim (Mq) 


Jin (Jkaw) 


Msec (Mo) 


R(r.) 


log(k2) 


0.30 








-0.872 


0.30 




0.30 


7 


-0.877 


0.30 


1 






-0.921 


0.30 


1 


0.30 


7 


-0.919 


0.30 


1 


0.30 


14 


-0.918 


0.30 


1 


0.30 


35 


-0.919 


0.30 


1 


0.30 


70 


-0.920 


0.30 


1 


0.15 


7 


-0.918 


0.30 


1 


0.60 


7 


-0.919 


0.30 




0.30 


3 


-0.868 


0.30 


1 


0.30 


3 


-0.916 


1.00 








-1.904 


1.00 




i.bo 


"7 


-1.918 


1.00 


1 






-1.924 


1.00 


1 


i.bo 


"7 


-1.923 


1.00 


1 


1.00 


14 


-1.926 


1.00 


1 


1.00 


35 


-1.935 


1.00 


1 


1.00 


70 


-1.930 


1.00 


1 


0.50 


7 


-1.932 


1.00 


1 


2.00 


7 


-1.939 


1.00 




1.00 


3 


-2.094 


1.00 


1 


1.00 


3 


-2.027 


1.00 


3 


1.00 


3 


-2.277 



''Ellipsis dots (. . . ) indicate missing values, because the 



model is either non-rotating or not binary. Rows with all 
numbers correspond to rotating binary models. 



than for those calculated with the standard models. This 
is due to the combined differences between the radius and 
the rotational inertia of each model. 

Concerning the effects on the internal structure of the 
stars, the models distorted only by tidal forces are in gen- 
eral similar to the standard ones and the models distorted 
only by rotation are similar to the rotating binary models. 
This is due to the fact that, for the adopted orbital sepa- 
ration, the rotation effects are, on average, more important 
than the tidal ones. Here, we concentrate in the differences 
between the standard models and the rotating binary mod- 



els, since this is the case of scientific interest for apsidal 
motion studies. 

In Fig. [2] we show the path followed by our standard 
models (solid lines) and our rotationally and tidally dis- 
torted models (dashed hues) for the masses of 0.09, 0.2, 0.3, 
0.4, 0.6, 0.8, 1.0, 1.2, 1.8, 2.0, and 3.8Mq. We can see that 
the distorted pre-MS evolutionary tracks have lower Toff 
than their standard counterparts, a s it is already kn own 
for stars distorted only by rotation. ISackmaiml (|1970( ) re- 
ported that the distorting effects of rotation affects late- 
type, main-sequence stars much more seriously than the 
early-type ones, and that this behavior could be approxi- 
mately understood by looking at the variation of the power 
law expressions of the nuclear-energy-generation rates along 
the main sequence. This trend is also observed when tidal 
and rotational distortions act together during the pre-MS. 

Fig. [31 in which we plot logA:2 against the stellar age for 
selected masses of our standard models and rotating binary 
models, shows the significant change in mass concentration 
during the pre-MS evolution, especially for higher masses. 
For ages less than 1 Myr the log k2 do not vary significantly, 
neither with time nor with mass (see Fig. [3]). For models 
with masses lower than 0.3 M©, that correspond to the mass 
range where stars arc almost completely convective, log k2 
remains roughly constant during the evolution. For masses 
greater than 0.3 M© the values of logfc2 are constant until 
a given age, after which they start to drop; the decrease 
of log ^2 starts earlier as the mass increases. It seems that 
log fc2 remains constant until the star develops a radiative 
core. 

The values of logfcj versus stellar mass for standard 
models and rotating binary models at the ZAMS are shown 
in Fig.m in which we see that the distorted models corre- 
sponding curves remain below those corresponding to the 
standard models. The values of log kj are more sensitive to 
the distorting effects for M<0.5Mq and for M>1.4M0. 
The values of log k2 , log fca , and log k^ remain roughly con- 
stant in the mass range 0.09-0.4 M©. As the mass increases 
from 0.4Mq to 1.5Mq, kj drop significantly (2-3 orders of 
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Fig. 2: Evolutionary tracks for standard and for rotating 
binary models. The model masses eorrespond to 0.09, 0.2 
0.3, 0.4, 0.6, 0.8, 1.0, 1.2, 1.8, 2.0, and 3.8 from bottom 
to top. 
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Fig. 3: The time evolution of log fc2 for selected masses of 
the standard and the rotating binary models. The gray line 
corresponds to the ZAMS. 



magnitude) and reach their minimum value at 1.5 Mq. For 
masses greater than 1.5 we note a parallel behavior of 
the kj. H87 pointed out that, for a given mass, there is a 
significant decrease in kj for increasing j , and that the rele- 
vance of including the higher order terms, when comparing 
with observations, can be judged from Eq. In our anal- 
ysis logfcs and logfc4 are really less important than logA:2 
in the mass range that he analyzed (0.5-32 Mq), but the 
same statement cannot be extended to less massive stars 
(M<0.5Mq): as can be seen from Table [2] and also from 
Fig. 21 is lower than k2 by more than one order of mag- 



Fig. 4: \ogkj as a function of log A/ for ZAMS models, for 
the the same set of masses as in Fig. [21 



nitude in the mass range 1.5-3.8 M©, but this difference 
decreases to less than one order of magnitude for stars less 
massive than 0.5 Mq. So, in the low-mass range, the as- 
sumption that the harmonics of order greater than j=2 
can be neglected, widely used when analyzing the apsidal 
motion of binary systems, seems not to be entirely justified, 
except maybe when computing the apsidal motion rate Co 
in which, besides a dependence on log kj , there is also a 
dependence on {ro/Ry^^ that will make the j ~ 2 term 
always dominate unless rp is very close to R. 

In Fig. m we show log/3 as a function of logM, at the 
ZAMS, for our standard and n on-standard models and als o 
for the non-rotating models bv lClaret fc Gimened (|l989b[ ). 
For masses abov e 1.4 Mq our model s predic t a lower value 
of log/3 than the lClaret fc Gimened (|l989bf ) ones; the val- 
ues predicted by our rotating binary models are even lower. 
Although these models have slightly different chemical com- 
positions, the comparison is still valid. A control made with 
a one solar mass mo del with the same initial ch emical com- 
position as those bv lClaret fc Gimened (|l989b| ) shows that 
the gyration radius at the ZAMS did not vary significantly. 

The time evolution of log (3 has the same behavior as the 
time evolution of log kj (Fig.jS]). Also, log /3 behaves roughly 
in the same way as logfcj as a function of logM (Fig. [4]), 
with a minimum at 1.5 Mq, corresponding to the change in 
the do minant energy source from the p-p chain to the CNO 
cycle (jClaret fc Gimened Il989b[ ) . This similarity is due to 
the known linear relationship between log kj and log /3, also 
seen in Fig. 1 by Jylotz (1952). 

The differences of the stellar radii as a function of stel- 
lar age and mass for the standard and distorted models are 
shown in Fig.|6l For a 0.5 Mq star, it can be seen that the 
stellar radii at the ZAMS of the rotating binary models 
are slightly smaller than those produced by the standard 
models. For the 1 Mq model, the situation is the opposite: 
the distorted models predict a slightly larger stellar radius 
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Fig. 6: Stellar radius and main axis of our standard models and rotating binary models. We report only the results for 
0.5, 1.0 and 2.0 M© models. The solid and dotted lines denote the radius of standard models and rotating binary models, 
respectively, while the dashed line stand for the radius along the main axis of the rotating binary models. 
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Fig. 5: Standard models (open circles) and rotating bi- 
nary models (crosses) in the log/3 - log M plane. The full 
triangles correspond to the results by IClaret fc Gimenea 
(|1989bD . 



at the ZAMS than the standard models; this difference be- 
comes still larger for the 2 M© m odel. The thresho ld mass 
for this transition is at ^O.TM©. ISackmaiml (|1970[ ) had al- 
ready noted this behavior, finding a threshold mass about 
1.3-1.5 M© and suggested that it could be related to the 
cross-over from the proton-proton chain to the CNO cycle 
that occurs around 1.5 M©, depending on the initial chem- 
ical composition. According to him, this behavior of the 
effects of rotation shows up in all physical quantities of a 
star. Fig.[n]also shows the stellar radius along the mean sys- 
tem axis, defined as the axis joining the two stars. For all 
models computed the mean axis is greater than the mean 
stellar radius during the whole evolution, the difference be- 



ing more evident for ages higher than ^lOOMyr (0.5M©), 
~20Myr (l.OM©) and -lOMyr (2.0M©). 

As a final remark, we verified that tidal effects act in 
the same way as rotational ones but on a smaller scale. 

6.2. Comparison with otiier works 

Computed values of internal structure constants for main- 
sequence stars are available in the literature from other 
authors. We checked the validity of our calculations by 
comparing our results of logfc2 at the ZAMS with those 
previously published ones. Our standard values of internal 
structure constants are in qualitative agreement with those 
by H87, CG89a and CG92, as well as our distorted logfc2 
relative to those given by C04, COS, C06 and C07. To bet- 
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Fig. 7: log k2 as a function of the stellar age for a 1 M© 
standard model, and for the initial chemical compositions 
(X, Z)= (0.7, 0.02), (0.07, 0.04), (0.7, 0.004), and (0.8, 0.02). 

ter compare our standard results with the previous ones. 
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we computed additional grids of 1 Mq standard models, 
with different values of the mixing length parameter (a~2.0 
and a=1.5) and four different initial chemical compositions. 
(X,Z)=(0.7,0.02), (0.7,0.04), (0.7,0.004), and (0.8,0.02)! 
From Fig. [3 where the time-dependence of log/c2 for a=1.5 
and the initial chemical compositions above mentioned are 
shown, one can have an insight on how the metallicity af- 
fects the value of the second-order internal structure con- 
stant: during the first 1 Myr the chemical composition does 
not alter logfc2, but it becomes important from this age on; 
the metal-poorer stars evolve to more centrally condensed 
configurations (lower values of logfc2) than their metal- 
richer counterparts. As it can also be seen from Fig. [71 using 
the initial chemical composition (X, Z)=(0.7, 0.04) (dotted 
lines) produces roughly the same effect on the evolution 
of logfc2 as using the initial chemistry (X, Z)=(0.8, 0.02) 
(long-dashed lines). 

In TableH we list the values at the ZAMS of the in- 
ternal structure constants obtained from our standard and 
distorted 1 Mq models and also from other works in the lit- 
erature. The values for our standard models were obtained 
with the same initial chemical compositions of Fig.[71 they 
can be compared with the values by, e.g., H87, CG89a and 
CG92. 

We start by noting that for the initial chemical compo- 
sition (X, Z)=(0.7, 0.02) and our IMq model pro- 
duces a lower value of logfc2 than that obtained by CG92, 
while for the same metallicity and q;=2.0, it produces a 
lower value of log k2 than those of CG89a and H87. Now, 
let us consider the a that fits the sun for each model, 
namely a=1.5 for this work and CG92, and a=2.0 for 
CG89a and H87. Still keeping the same metallicity, our 
models produce the most centrally condensed 1 Mq star, 
followed by that of H87, CG89a and CG92, in this or- 
der. These differences between our values of log fc2 and 
those by CG89a and CG92 can be attributed to the use 
of different opacities, as they used Los Alamo s Opacity 
Libraries and OPAL ([Rogers fc Iglesiaj . Il992f ) opacities 
respe ctively, while we used m ore recent OPAL opacities 
from llglesias fc Roger j (|1993[ ): in addition, CG92 consid- 
ered mass loss in their models. It is worth noting that al- 
though CG92 used more updated opacities than CG89a, for 
the case of (X, Z)=(0.7,0.02) the former authors obtained 
a value of ^2 for their 1 M© model at the ZAMS 1.34 times 
greater than that obtained by the latter ones. On the other 
hand, the CG92 models produce lower values of fc2 than 
those of CG89a for masses greater than 1.1 M©. 

For the remaining initial chemical compositions, we can 
only compare our 1 M© standard model with those by 
II87. For the case of Q!=2.0, we notice that H87 obtained 
less mass-concentrated models for (X, Z) = (0.7, 0.004) and 
(X,Z)=(0.8,0.02) than us, while for (X, Z)=(0.7, 0.04) it is 
the opposite. We can attribute this fact to two reasons: (1) 
although H87 obtained the lowest values of k2 for one of 
the four meta licities used in this com parison, he used older 
opacities from lCox fc Stewart! (jl969t ): and (2) our best fit to 
the sun is obtained with a=1.5, not with Q!=2.0. Regarding 
this last point, if we turn our attention to the models for 
which the value of a can reproduce the solar radius at the 
solar age for this work and a=2.0 for H87), we see 

that our models again produce a more centrally condensed 
star than those of II87, independent of the chemistry. 

Finally, we compare the log(A:2) values from our dis- 
torted models with those of C04, C05, C06b and C07. 



Table 4: Values of log k2 for a 1 M© star model at the 
ZAMS, obtained with our models (standard and distorted) 
and those previously published in the literature. 



(X,Z) 


Qmlt 


Reference 


log(fc2) 




2.0 


H87 


— 1.768 








1 747 


(0.7, 0.02) 


2.0 
1.5 


(standard) 
CG92 


-1.804 
-1.619 




1.5 


(standard) 


-1.975 




2.0 


H87 


-1.614 


(0 7 041 


2.0 
1.5 


(standard) 
(standard) 


— 1.548 
-1.628 




2.0 


H87 


-2.035 


(0.7, 0.004) 


2.0 
1.5 


(standard) 
(standard) 


-2.181 
-2.281 




2.0 


H87 


-1.578 


fO 8 02) 


2.0 
1.5 


(standard) 
(standard) 


— 1.688 
-1.776 


(0.7, 0.02) 


1.68 
1.68 
1.68 


C04 
(standard) 
(distorted) 


-1.583 
-1.859 
-1.895 


(0.748, 0.004) 


1.68 
1.68 
1.68 


C05 
(standard) 
(distorted) 


-1.840 
-2.112 
-2.145 


(0.73, 0.010) 


1.68 
1.68 
1.68 


C06b 
(standard) 
(distorted) 


-1.682 
-1.981 
-2.004 


(0.64, 0.04) 


1.68 
1.68 
1.68 


C07 
(standard) 
(distorted) 


-1.523 
-1.852 
-1.866 



In order to perform suitable comparisons, we computed 
additional 1 M© models with the same physical inputs 
used in those works, namely q;=1.68 and initial chemical 
compositions of (X, Z)=(0.7, 0.02), (X, Z)=(0.748, 0.004), 
(X,Z)=(0.730, 0.010), and (X, Z)=(0.64, 0.04); the corre- 
sponding values of logfc2 at the ZAMS are shown in the 
lower part of Table [H We note that although the effects 
of tides and rotation are a lso considere d in the models by 
Claret and collaborators (jClaretl [l99i: Claret 2008, per- 
sonal communication), their values presented in Table |3] 
correspond to isolated, non-rotating stars; for comparison 
purposes, we also present in it our values calculated for iso- 
lated, non-rotating stars with the same masses and chemi- 
cal composition as well as the corresponding ones including 
tidal and rotational distortions. For all chemical composi- 
tions, our models again produce more mass-concentrated 
stars than those by C04, C05, C06b, and C07. For the 
mass range of this comparison, these differences in the 
logfc2 values can be due to the presence of distortion ef- 
fects in our models and also to th e use of a different e qua- 
tion of state: while we use OPAL ()Rogers et al.l . Il996f) and 
iMihalas et all (Il988l) EOS, those autho rs use the CEFF 
(jChristcnsc n-Dalsgaard fc Dappenl . Il992[) EOS. 

Fig. [5] presents a graphical comparison between the val- 
ues of log k2 obtained in this work and those obtained by 
H87, CG89a, CG92, C04, COS, C06b, and C07. This figure 
gives an easier overview of TableSl 

7. Comparison between theory and observations 

Double-lined eclipsing binary systems are good candidates 
to test evolutionary models, but very few systems are as ad- 
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Fig. 8: A schematic comparison between logfc2 values ob- 
tained with our standard models (Sq;=1.5 and Sa=:2.0); 
our distorted models (Da=1.68) and those obtained by 
other authors. The values plotted are the same as in 
Tabled! and the corresponding initial chemical com- 
positions are the following, labeled A to G respec- 
tively: (X,Z)=(0.7,0.004), (0.7,0.02), (0.7,0.04), (0.8,0.02), 
(0.748,0.004), (0.73,0.01), and (0.64,0.04). 



equate as EK Cep (P=4'*.42). The mass and radius of its pri- 
mary component are Mi=2.029±0.023Mo and Ri = 1.579± 
O.OO7R0, while, for the s econdary, M2= l. 124 ± 0.012 Mf;, 
and R2=1.315 ± 0.006 Rw (|Claretl . l2006ah . As mentioned by 
IClaretl (|2006a[ ). EK Cep has accurate determination of ab- 
solute parameters, at least concerning masses and radii (see 
Tablets]); its secondary component is a pre-MS star, the ap- 
sidal motion presented by the system has a high relativistic 
contribution, the less massive component has its Lithium 
abundance determined observationally, the metallicity of 
the binary is evaluated, and, further, the rotational veloc- 
ity of each component is observationally measured. Besides 
that, EK Cep is the only known pre-MS system with mea- 
sured apsidal motion. For all these reasons, we choose EK 
Cep to test our new models. 



Table 5 : Abso lute dimensions of EK Cep (jClaretl . l2006at 
lTomkinl . fl983l) . 



Parameter 



Mass (Mq) 
Radius (R©) 
log(5') (cgs) 
log (L/Lo) 
log {T.b){K) 
Vrot (km/s) 



Primary 



Secondary 



2.029 ± 0.023 
1.579 ± 0.007 
4.349 ± 0.010 
1.17 ± 0.04 
3.954 ± 0.010 
23±2 



1.124 ± 0.012 
1.315 ± 0.006 
4.251 ± 0.006 
0.19 ± 0.07 
3.756 ± 0.015 
10.5±2 



EK Cep w as discovered as an eclipsing binary 

by IStrohmeierl ()1959D from photographic observations. 
Photometric elements and a spectroscopic o rbit of the sys- 
tem were first presented by lEbbighausenl (|l966al lbl). and 



further revised by other authors. The apsidal motion of EK 
Cep was first reported in lKhaliuUinI ([l983i) . [^mki3 ^9M) 
determined the masses and the radii for its primary and 
secondary components. He also noted that the secondary 
is oversized in comparison with a main-sequence star with 
the same mass and supposed that it might still be con- 
tractin g towards the main sequence. iHiU fc Ebbighausenl 
concluded that bo th components are zero- age main- 



sequence obje cts, despit e TonikinI (Il983l ) suggestion about 
the secondary. iGimenej (|1985[ ) obtained a good agreement 
between theoretic al and o bserv ationally determined apsi- 
dal motion rates. iPoppeil (|l987t ) reported some anomalies 
in the secondary of EK Cep, such as low effective gravity 
and temperature as well as the excess radiation in the blue 
band, that appear to be consistent with the hypothesis of 
the pre-MS nature of this star. Fr om high-resolution spec- 
troscopy in the Lil A6708 A region. iMartin fc Rebolol (|1993[ ) 
determined the lithium abundance of EK Cep B and pro- 
vided new eyidence that it has not settled onto the ZAMS. 
IClaret et al.l (|1995H ) compared the observed parameters of 
EK Cep with theoretically predicted values; they derived a 
common age for the system around 2x 10^ yr and confirmed 
the fainter component as a pre-MS star, while the more 
massive companion is in the beginning of the Hydrogen- 
burning phase. They estimated a Newtonian apsidal motion 
rate that is in agreement with the observations, considering 
a predicted relativistic contribution of about 40%. Those 
authors also found that the lithium depletion computed by 
th eir models is consistent w ith the abundances determined 
bv lMartm fc Rebolol (|1993D . 

The evolut ionary statu s of EK Cep was studied by 
other authors. lYildid ()2003f ) modeled the component stars 
by invoking a rapid ly rotating core for the primary. 
iMargues et al.l (j2004f ) investigated the role of overshoot- 
ing on t he mo deling of pre-MS evolution of the secondary. 
IClaretl (|2006aD presented the most recent analysis about 
EK Cep; due to problems with the empirical determination 
of the effective temperatures of the component stars, he 
adopted the effective temperature ratio (TR=Toff,2/Teff_i), 
which is better determined from the light curve analysis 
than their absolute values. Inconsistencies found in the pho- 
tometri c distances for both components support this ap- 
proach. [CWi3 (j2006af ) used other constraints in his analysis 
as the radii, apsidal motion and lithium depletion. With a 
rotating model (assuming local conservation of the angular 
momentum), a=1.4 and (X, Z)=(0.7075, 0.0175), he fitted 
the radii and TR in the same isochrone (24.2x10^ yr). 

As emphasized by IClaretl (|2006al ). before computing 
stellar models for a given star it is fundamental to define 
clearly which observational parameters will be used as con- 
straints. Following the same approach as lClaretl (|2006a[) . we 
adopt the masses, radii and TR, leaving the apsidal motion, 
rotational velocities and lithium depletion to be used as ad- 
ditional constraints after obtaining an acceptable solution. 
We then computed rotating binary models with metallici- 
ties as close as po ssible to the solar metall icity, following 
the conclusion by iMartm fc Rebolol (|1993D that the EK 
Cep secondary has a metal contents typical of an young 
disk, solar-type star. Using an initial chemical composition 
of (X, Z)=(0. 67, 0.017) and the mixing length parameter 
Q!=1.5, we obtained a model that reproduces the radii and 
TR for EK Cep within the uncertainties, as shown in Figs.jSj 
and llOi deriving two possible age intervals for the system, 
15.5-16.7x10^ and 18.9-19.3x 10^ yr. The input values we 
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Fig. 9: Radii predicted by our binary rotating models for 
EK Ccp components determined masses. The additional 
tracks refer to the maximum and minimum mass of each 
component, according to their errors. Full lines denote the 
primary while dotted ones denote the secondary. Note that 
there is an acceptable agreement between radii and effec- 
tive temperature ratio (see Fig. fTO]) for two age intervals 
(vertical lines). Horizontal lines represent the error bars in 
radii determinations. 
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Fig. 10: Effective temperature ratio (TR=Toff,2/Toff,i) pre- 
dicted by our binary rotating models for EK Ccp com- 
ponents determined mass. The additional tracks refer to 
the maximum and minimum mass of each component star, 
according to the errors in their determination. The verti- 
cal lines correspond to the same age intervals as in Fig[9l 
Horizontal lines represent the error bars in TR determina- 
tions. 



used are different from those by IClaretl HqOGJ) since we 
are using different stellar models, with different physics. By 
using his inputs, i.e. (X, Z)=(0.07075, 0.0175) and a = 1.4, 
we were not able to reproduce, simultaneously, both stellar 
radii at the same age interval; this is due to the fact that 
each model needs slightly different inputs for adjusting to 
a same set of observed data. 

According to our radii and TR analysis, the two derived 
age intervals allow a reasonable adjustment of the observa- 
tional data; we used additional data (rotational velocities of 
each component, lithium depletion and apsidal motion con- 
stants) to constrain even more the system's age estimate. 

We first tried to reproduce the observed rotational ve- 
locites of EK Cep primary and secondary components at at 
least one of the derived age intervals; in order to do so, we 
needed to use an initial angular momentum for each com- 
ponent lower than that prescribed by Eq.[72l In Fig.fTT] we 
show the time evolution of the rotational velocities for EK 
Cep determined masses; the vertical lines correspond to the 
two age intervals that fit simultaneously the stellar radii 
and temperature ratio, while the horizontal ones delimit 
the observed EK Cep rotational velocities. From that fig- 
ure, one can see that the older age interval does not match 
the observed rotational velocity of the secondary star. 

The lithium depletion of our rotating binary models 
is shown in Fig.[T21 in which the Li contents is plotted 
against the stellar age. We starte d from an initial lithium 
abun dance of log (Li/H)=3.1 (as lD'Antona fc Montalbii] . 
l2003f ). As expected, we do not find any depletion for the pri- 
mary during its pre-MS evolution due to its higher mass; 
for the secondary, we find a depletion of about 0.28 dex in 
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Fig. 11: Time evolution of rotational velocities along the 
pre-MS for EK Cep determined masses. Horizontal lines 
represent the error bars in rotational velocity determina- 
tions. Curve triplets and vertical lines defined as in Fig. [5) 



both possible age intervals for EK Cep. This corresponds 
to a lithium abundance of log (Li/H)=2.82±0.03, consistent 
with the surface value o f log (Li/H)=3.1±0.3 measured by 
iMartm fc Rebold (|1993[ ). The lithium depletion does not 
allow us to choose between the two possible age intervals 
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Fig. 12: Time evolution of lithium abundances for EK Cep 
determined masses according to our rotating binary models. 
Curve triplets and vertical lines defined as in Fig.[9l Note 
that the curves corresponding to the primary mass and to 
its error bar collapse to a single horizontal line. 
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Fig. 13: The second-order apsidal motion constant as a func- 
tion of the stellar age for EK Cep, as obtained from our 
rotating binary models. Curve triplets and vertical lines 
defined as in Fig.l^] 



for the system, since it does not vary in the whole interval 
from 10 Myr to 50 Myr. 

Finally, we investigated the apsidal motion rate of 
EK Cep with the same rotating binary models. The val- 
ues for the orbital eccentricity (0.109 ± 0.003 degrees) and 
the orbital inclination (89.3 ± 0.1 degrees) were taken from 
iPetrova fc Orlovl (jl999[ ). The anomalisti c period of EK Cep 
is Pan = 4.4278062 ± 0.0000005 days (|Claretl . [2006ah . In 



Table 6: Summary of the quantities related to the apsidal 
motion. 

log ^2 U (years) uj (o/cyclc) 



Our models (All) 
Our models (AI2) 
Observed 
Claret (2006a) 



-2.1 



m 



0.00j)|| 
0.000|^ 
0.00^5^ 
0.00095 



Fig. [13] we show the variation of log(A:2) for the observed 
masses (with their errors) as a function of age. As before, 
the vertical lines indicate the age intervals that fit both 
the radii and effective temperature ratio of EK Cep com- 
ponents, and from which we took the values of the internal 
structure constants to derive the apsidal motion rate. 

For the first age interval, 15.5-16.7Myr (All), our mod- 
els result in ks values of 0.00261±0.00001 and 0.007±0.003 
for the primary and the secondary components, respec- 
tively. From Eq. (fTO)) we obtain log fc2thco=— 2.3±0.2 as the 
mean value of the second-order apsidal motion constant 
of the system. The Newtonian contribution of the apsi- 
dal motion rate is obtained from Eqs. ([3]) and p^ . re- 
sulting in wjv=(0.00033±0.00011) "/cycle. We also calcu- 
lated the relativistic contribution of the advance of peri- 
astron, a;i^=(0.000434±0. 000005) "/cycle, which, based on 
our models, corresponds to about 49-68% of the total rate. 
From these values, our models result in an apsidal motion 
rate of cj=(0.00076±0. 00012) "/cycle and, consequently, an 
apsidal period of Uthco=(5800±800) yr; these values are 
comparable with the observed ones. However, for the sec- 
ond age interval, 18.9-19.3 Myr (AI2), we were not able to 
reproduce the observed apsidal motion of EK Cep by using 
our rotating binary models. As can be seen in TablelHl the 
apsidal motion quantities we found in AI2 (log/c2, apsidal 
period U and apsidal motion rate uj) are not in agreement 
wit h the o b served values, even considering the error bars. 

Claret (l20()6a ) used results of a series of works 

(iKhafiuUinl Il983bl : ICimenez fc Margrave! 119851 : 



iHiU fc EbbighausenI [T984I: IClaret et all ' Il995bt and 

others) to derive the observed apsidal motion rate of 
this system. He reported the observed _mean value of 
the internal structure constant as log fc2=— 2.09±0.09, 
which is equivalent to an observed apsidal motion of 
about tjobs=(0.00097±0.00015) "/cycle, after applying the 
relativistic correction. This apsidal motion rate produces 
an apsidal period of Uobs= (4500±700) ye ars. From Table[6l 
where our results, those bv lClarej (|2006a[ ) and the observed 
determinations are summarized, it can be seen that our 
predicted values for All are in g ood agreement with the 
observed ones (as well as those by IClaretl . l2006a| ) , and the 
differences lie within the errors. 

At first glance, we could be tempted to conclude that 
the second age interval AI2 gives a better agreement for 
the EK Cep system than the first one due to the following 
reasons: (1) at AI2, the radius of the secondary component 
constrains the age of the system to a much narrower age in- 
terval, as shown in Fig.[9l and (2) the TR curves produced 
by our models lie very well inside the AI2 age interval, but 
are fitted just marginally for the earlier All age interval 
fFig. fTO)) . However, in that age interval our rotating binary 
models were not able to reproduce neither the rotational 
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velocity of the secondary star fFig. [TT|) nor the apsidal mo- 
tion rate for the system (Tablc[S]), while the opposite did 
happen for All. Therefore, All seems to be more suitable 
than AI2 as the age interval for investigating the evolution- 
ary status of EK Cep. For these reasons, with regard to our 
models, we suggest that the age of EK Cep is between 15.5 
and 16.7 Myr. 



8. Conclusions 

We computed stellar evolutionary models that take into ac- 
count the combined effects of rotation and of tidal forces 
due to a companion star, obtaining values of the internal 
structure constants for low-mass, pre-MS stars. Our ap- 
proximation for binary model calculations (with effects of 
tidal forces) do not include the pseudo-potential centrifu- 
gal terms arising from the orbital motion of the components 
around the center of mass of the system. For all sets of mod- 
els, namely standard, binary, rotating and rotating binary 
models, we tabulated the internal structure constants and 
the gyration radii for ZAMS models and their time evo- 
lution. Distorted models result in more mass-concentrated 
stars and produce larger gyration radii at the ZAMS than 
standard ones. 

The non-standard evolutionary tracks are cooler than 
their standard counterparts, mainly for low-mass stars. 
Regarding the internal structure of our stellar models, we 
verified that tidal effects act in the same way as rotational 
ones but in a smaller scale; the relative importance of these 
effects on the apsidal motion constants depend mainly on 
the orbital separation and the star's rotation rate. Besides, 
rotationally and tidally distorted models produce some ef- 
fects in the physical quantities of a star, in comparison with 
the standard models, that are in one direction, for masses 
below ^0.7 Mq, and in the o pposite dir e ction, for masses 
above this threshold. Though TSackmaniil |l970) associated 
this behavior in distortion effects with the transition from 
the p-p chain to the CNO cycle, this transition occurs in 
a mass range about 1.3-1.5 Mq, considerably above the 
threshold we observed with our models for the change of 
behavior of the physical quantities of the stars. 

We also found that, for masses lower than 0.5 Mq, the 
relative importance of the second-order internal structure 
constants over those of higher order is much lower than 
for the mass interval of 1.5-3.8M0; hence, in the low-mass 
range, the usual assumption that the harmonics of order 
greater than j=2 can be neglected seems not to be entirely 
justified except maybe for computing the apsidal motion 
rate uj. 

Our results on internal structure constants were com- 
pared with those available in the literature and found to 
be compatible with them. The k2 values obtained from our 
standard models are smaller than those last published by 
H87, CG89a and CG92, except for a given model by H87 
with the inputs (X, Z)=(0.7, 0.04) and a=2.0. Our rotating 
binary models produce internal structure constants even 
smaller, resulting in more mass-concentrated configurations 
than the models by C04, C05, C06, and C07. These com- 
parisons were made for representative 1 Mq models, with 
the same initial chemical compositions of those works. 

Using our set of evolutionary tracks for rotating binary 
models, we also investigated the evolutionary status of the 
interesting double-lined eclipsing binary system EK Cep. 



Its primary, a 2.029 M© star, seems to be in the hydrogen- 
burning phase, and its secondary, a solar-like star (1.124 
Mq), is confirmed as a pre-MS star. By using a model with 
an initial chemical composition of (X, Z)=(0. 67, 0.017) and 
a mixing length parameter of a=1.5, we reproduced stellar 
radii and the effective temperature ratio of EK Cep in two 
different age intervals. We also followed the lithium con- 
tents during the pre-MS evolution of both components; as 
expected, we do not find any significant depletion for the 
primary, while the Li depletion for the secondary agrees 
with the observed values within the uncertainties for both 
age intervals. However, for the later age interval of 18.9- 
19.3 Myr, we were not able to fit neither the observed rota- 
tional velocity of both components nor the observed apsidal 
motion rate of the system. On the other hand, all observed 
quantities used in our analysis were reproduced in the ear- 
lier age interval of 15.5-16.7 Myr. Although being a broader 
range for the age of EK Cep and introducing greater error 
bars to the apsidal motion quantities, it seems to give a 
more confident age for the system. 

In this paper, we presented a first attempt to intro- 
duce the combined structural effects of tides and rotation 
in the ATDN evolutionary code. Even with simple assump- 
tions about the secondary component of a binary system, 
we could verify the importance of those combined effects on 
stellar structure and evolution. Work is in progress, where 
we remove the assumption that the disturbing star can be 
treated as a point mass, include the pseudo-potential cen- 
trifugal terms caused by the orbital motion and consider the 
simultaneous evolution of the two components of a binary 
system. 
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